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1. INTRODUCTION 


I,] Fundamental Philosophy of Physiological System Modelling 

The importance of physiological system modellihg within the area of biomedical research 
is becoming more evident. In particular, the application of control theory which involves 
the mathematical description of element functionings and interrelationships between elements 
has made possible the simulation of very complex physiological feedback control systems. 

The mathematical relationships range from algebraic formulations and ordinary linear 
differential equations to nonlinear differential -difference equations and stochastic processes. 

Generally, physiological systems are very complex and comprise many interrelated 
variables. In the initial phases of model formulation verbal descriptions are used to formulate 
the expressions that are based upon experimental physiological data or assumptions. This 
approach of verbally explaining the functioning of the physiologi cab system becomes very 
unwieldy as the system's complexity increases. For this reason specific variables and symbols 
are defined and related to each other by mathematical expressions. Once the system is 
defined in logical mathematical terms the analysis process is then adapted to various computa- 
tional devices, i.e. digital, analog, and hybrid computers. By altering parameter and 
functional relationships the model's responses are correlated with experimented data. At 
this stage in the model's development, the refinements made by the investigator determines 
how well the model faithfully describes the actual system. 

According to Yamamoto and Raub (1) the best means of categorizing the present stage of 
development of respiratory control system models would be to desciibe them as being in a 
transition stage from verbal to mathematical statements. There are still many physiologically 
related lespiratory functions that are not easily desctibed mathematically. One necessary 
component of the oveiall plan to rectify this situation is the establishment of experimental 
efforts that provides quantitative measures instead of qualitative measures. From this point 
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of view perhaps the respiratory system and the related simulation efforts are even more firmly 
established than for other physiological systems. In particular, the control mechanism 
associated with respiration and its interaction with the mathematical description of the 
controlled component is defined by several sources (2-8). 

Prior to deciding upon a particular physiological system model one should utilize the 
expertise found in the literature and correlate the features of the published models with the 
goals of the immediate project. This is precisely the approach that was established in this 
research effort involving the respiratory control systems. 

1 .2 Various Respiratory System Models 

The system modelling effort that is discussed in the following report is one related to the 
respiratory control system. Several models have been developed for the respiratory control 
system emphasizing particular features which make them adaptable for specific experimental 
and simulated environmental conditions. 

A few brief comments about some of the respiratory modelling efforts and their salient 
characteristics follow. Alveolar ventilation responses were of prime importance for the 
output responses of a model developed by Belivilfe et al. (9). Frequency, amplitude, and 
mean level of inhaled CO 2 were altered in a manner corresponding to sinusoidally varying 
CC >2 inhalation allowing a frequency resrxsnse analysis of the respiratory system. 

The phenomenon of pulmonary mechanics as related to various forms of inspired-expired 
flow systems has been researched by Fry (10). Here variables related to physiological 
and pathological features of the air passageways are better defined. Elastic behavior of the 
lungs, airway compliances, and resistances weie consideied by Mead (11) in a model possessing 
characteristics similar to those of Fry's model. 
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Sai del et al. (12) developed a lumped parameter simulation of the pulmonary gas 
transport system with the ultimate goal of analyzing parameter changes. The model included 
five variable-volume compartments approximating the airway and alveolar regions and a 
constant-volume compartment representing the capillary bed. Extensive efforts were taken 
to mathematically describe the mass balance of the gas compartments and blood compartments. 

Breathing patterns related to CO 2 oscillatory exchanges are presented in a model 
developed by Yamamoto and Hori (13). A linear controller equation is used to describe, 
indirectly, the dependence of tidal volume and rate of breathing upon cerebral PcC>2' 

Exercise levels (workload) are very important inputs for the proposed research. Analysis 
of maximum ventilation related to chronic airway obstruction and exercise has been conducted 
by Pierce et al. (14). This particular research was not intended to be a respiratory control 
system model study; however, the results of the endeavor are very applicable fo modelling 
of airway passages and to the investigation of the respiratory-exercise relationships. 

A more detailed discussion of previous efforts involving exercise-respiratory phenomena 
5s given in Section 1.3. Both humoral and neural control mechanisms are presented with 
supportive experimental work by several investigators (15-23). 

Milhorn and Brown's steady-state human respiratory system (5) is very closely related 
to Grodins 1 respiratory control system (2). In this particular model by Milhorn, a controlling 
system and a controlled system form the basis from which mass balance equations are written 
describing the exchange mechanisms for CO 2 and 02- The controlling equation of Milhom's 
model is an important fundamental component and is similar to one of the forms in Grodins 1 
mode! except for the void in exercise dependence and the monitoring of CC>2 and 02 instead 
of compartmentaf H + concentrations. 
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Closely associated with Milhorn and Brown's model is a study of pulmonary capillary 
gas exchange by Milhorn and Pulley (24). This work portrays a more detailed version of the 
capillary gas exchange by introducing anatomical parameters into the model. 

Assuming a second order representation for the mechanical portion of the respiratory 
system Hiiberman et al. (25) has utilized the Fourier series analysis and phasor notation 
to deteirmne the modulus of impedance, phase angle, compliance, and resistance of the 
respircfoty system. In doing so this research has developed a technique for adjusting 
parameters associated with respiratory mechanics. Although this research effort is not 
directly related to the present research the evaluation of parameters for specific environmental 
conditions might be useful in extending the proposed model's capabilities. 

As outlined in the following sections of this report, Grodins 1 respiratory control system 
model is the one to which modifications are proposed. Since the pH of each compartment of 
fhis proposed model is an important monitored variable the research by Hermansen et al . (26) 
relating exercise and blood pH could prove to be a valuable piece of supportive research 

S 

in determining a level of maximal exercise for experimental subjects and in turn for the 
simulation efforts. 

A comparison of several respiratory system models' characteristics and capabilities are 
presented by Yamamoto and Raub (1). This particular review emphasizes that all models 
possess unique features but are categorized into three groups (a) Models related to those of 
Gordins, (b) Models related to Horgan's efforts, and (c) Models related to Milhorn's 
research. Aiso this review cites several other models giving their comparisons to the above 
three groups. 
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1.3 Theoretical Aspects of Respiration and Exercise Interrelationships 

It is a well established fact that exercise plays a very important role in the control of 
respiration. As stated by Strippoli ( 27 ) three significant variables that must be considered 
in determining the degree of respiratory system dependency upon exercise are: 

(A) Intensi ty of exercise, 

(B) Type of exercise, and 

(C) Environmental conditions 

The following material which explains the theories of respiratory control involving 
chemical, neural, and neural-chemical interactions was reported by Strippoli ( 27 ). The 
first approach utilized in explaining exercise-respiratory control interactions involved a 
simple explanation of decreasing arterial C>2 and increasing arterial CO2. The chemical 
reaction is defined as 

C 0 2 + H 2 0^£H 2 C03^H + + HCO3 - . (1 - 1) 

At this point uncertainty prevailed as to whether the response was caused by CC>2 itself 
or the H + concentration released by the hydrated form of CO2, carbonic acid (H2CO3). 
However, experimental evidence indicated that neither CO2, anoxia, or the combination 
of high CO2 and low C*2 levels would produce the magnitude of the ventilation response 
observed during exercise. This fact stimulated research efforts to justify yet another form of 
stimuli during exercise. The humoral, neurogenic, and neurohumoral theories are discussed 
in the following paragraphs. 

(A). Humoral Theories 

Three chemical substances are known which are related to exercise by perpnea . These 
are the arterial gas tensions of CO2, O2, and the H 1 concentiation . There is also a possibility 
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of a release of chemical substances formed in the muscles which affect the respiratory response, 
but these are not well documented. A brief discussion of the effects of the three forenamed 
chemical substances are included here. 

1. Carbon dioxide tension, 

The CC*2 production is increased during muscular exercise. This increase should 
lead to an increased concentration of CC>2 in the arterial system sufficient to stimulate 
the respiratory center so that the response (hyperventilation) would lead to a decreased 
blood CO 2 level. This regulation of respiration process would be similar to the 
operation of a feedback system. However, the blood CO 2 level increases very little 
during exercise, certainly not of sufficient magnitude to explain the observed change in 
ventilation. 

2. Oxygen tension, PQ 2 

One would expect to see a slight decrease in blood P 02 during exercise. However, 
this decrease would not amount to more than a few mm Hg, since the saturation of 
arterial blood with O 2 is not greatly affected by exercise. Oxygen tension of approximately 
60 mm Hg would be required 'o produce a realizable stimulus from the chemoreceptors 
in the aortic arch and carotid bodies (15). 

3. H + concentration 

An increase in H + concentration of the arterial blood due to CO 2 accumulation 
has never been observed in mild exercise. Only during severe exercise does the H^ 
concentration of the arterial blood significantly increase. This increase is caused 
by large amounts of lactic acid produced by the exercising muscles during severe exeicise 
when anaetobic metabolism accounts for the majority of eneigy. 
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Even if all of the humoral mechanisms described above were functioning simul- 
taneously the increase In ventilation (l/min) observed during exercise could not be 
justified. Within the limits of aerobic metabolism, even if there are no appreciable 
alterations in blood chemistry, there is a linear reldtionship between steady-state 
ventilation and oxygen uptake. 

(B). Neurogenic Theoiies 

Krogh and Lindhard (16) were the first to interpret the sudden increase of 
ventilation at the beginning of exercise on a neural basis. They postulated that impulses 
from the motor cortex of the brain were sent to the respiratory center to effect a change in 
ventilation. In a different study, experimentation illustrated that the passive movements 
of the limbs elicited a ventilatory response (1). It was then thought that impulses from 
mechanoreceptors in the muscles, tendons and joints transmitted to the respiratoiy center 
were the cause of the response. Comroe and Schmidt (17) found that the mechanoreceptors 
involved were limited only to the joints. 

While performing cross circulation experiments on dogs, Kao (18) showed that 
the hyperventilation accompanying muscular activity is limited to the exercising neural 
dog. The humoral dog whose head was perfused by blood coming from the active limbs of 
the neural dog showed no change in ventilation. Thus, the increased ventilation of the 
neural dog was not due to the chemical stimulus of the respiiatory center, since there were 
no changes in the blood of the humoral dog. 

From these experiments and the research of Jensen et al . (19) it is established that 
at the onset of exercise, pulmonary ventilation increases abruptly in a much shorter time 
period than the circulation time; that is, before chemical substances formed by the exercising 
muscles can reach the respiratory center. 
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Neurogenic controls olone are inadequate to maintain the delicate chemical 
balance of the blood during exercise at the same level as when the body is at rest. Therefore, 
both.humoral and neurogenic mechanisms are involved for the most effective control of 
respiration during muscular exercise. 

(C). Neurohumoral Theory 

Dejours (20) lecorded ventilation and arterial CC >2 tension, P q £-q , on ° s ' n 9^ e 
breath basis to detect the very fast changes in ventilation at the initial and terminal phases 
of exercise. The ventilation increases abruptly at the beginning of exercise, causing a fall 
in P a (X>2* The ventilation then increases more slowly and reaches a steady-state level in 
a few minutes. Simultaneously, the increases and returns to the pre-exercise level. 

Corresponding changes take place at the end of exercise. A sudden decrease in ventilation 
causes an increase in Pq^q^ followed by a slow decrease in P a CC>2 re sting value. 

Since the neurohumoral theory seems to be adequately justified by experimental 


results, it has been utilized in the description of exercise hyperpnea. A representation 
of ventilation using the neurohumoral approach is shown in Figure 1. 

The modification of Grodins 1 model which simulates the exercise phenomenon 
utilizes a form of both neural and humoral control (28) with Figure 2 illustrating the stimuli 
paths. However, there are still some uncertainties as to the degree of neural contra! in the 
on- and off-transient responses in ventilation (21). 

1.4 Comments on Physiological Systems' Interactions 


As the project progresses there must be a concerted effort in establishing the criteria 
for the interaction between the four physiological systems; respiratory, cardiovascular , 
thermal, and the renai/endocrine/body fluid systems. Research such as that ol Albeigoni 
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FIGURE 1 . Time course of ventilation during exercise 
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et al . (29) which attempts to define the significant variables used in the interfacing model 
of the circulatory and respiratory systems and the recently published research effort of 
Guyton et al. (30) will help in the formulation of these criteria. Guyton's research may 
prove most valuable since it attempts to integrate sofne of the dynamical interactions of 
various physiological subsystems; i.e. circulatory dynamics, pulmonary dynamics, regional 
blood Flow control and its relationship to P 02 ' tissue fluids and pressures, heart rate, and 
stroke volume. 

To fully appreciate the transition between model choices or the potential capabilities 
of a model an understanding of the preceding references is mandatory. This is certainly 
not an exhaustive listing, but several references have excellent bibliographies which makes 
it an adequate state-of-the-art listing for this research. 
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2. TERMINOLOGY 

There are many terms and symbols used in the respiratory control system simulations 
of Grodins and Milhorn which are common to both. The terms involved with these 
simulations are related to each other and in turn compared to orher simulations in a summary 
by Yamamoto and Raub (1). .However, there are some specialized terms which ha/e 
particular meanings when used only with their lespective models. 

Throughout the following report, including the study reports (See Appendices 6.1 and 
6.2), new terms are defined when they are used. A listing of the more common terms are 
presented in Tables 1-3 of this section. 

2.1 Terminology Associated with Grodins' Models 

Terms associated with the basic Grodins 1 respiratory control system are presented in 
Table 1 (2). 



Table 1 . Symbols and Terms Used in the Basic Grodins 1 
Respiratory Control System Simulation 

Symbol * 

Description 

Units 

a i 

Solubility coefficient for gas 
(j =C 02 , C> 2 , N 2 ) in blood 

i iters/I iter 
blood/atm, 37C 

a;j 

Solubility coefficient for gas 

(i = CC> 2 / C> 2 , N 2 ) in 
brain (i = B), 

tissue (i =T), cerebrospinal 
fluid (i = CSF). 

liters/liter 
blood/atm, 37C 

B 

Barometric pressure 

mm Hg 

(BHC03)j 

Bicarbonate content of blood 
(i = b), tissue (i ~T), cerebrospinal 
fluid (i = CSF), brain (i - B) . 

liteis CO 2 / 
liter i, 37C 


A Note: The symbol Ajj is defined as follows. A designates variable class, i designates model 
location of the variable, j designates The chemical species of the vaiiable. 
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Symbol 





(Hb) 

k 


Description 


Concentration of j (gas, Hb02, H 1 ) 
in blood at lung exit (i - a), 
in blood at brain entrance (i -aB), 
in blood at tissue entrance (i ~ at), 
in brain (i - B), 

in cerebrospinal fluid (i CSF), 
in tissue (i — T), 

in blood at lung entrance (i - v), 
in blood at carotid body (I ■= ao), 
in blood at brain exit (t = vB), 
irt blood at tissue exit (i ~ vT). 

Diffusion coefficient of 

j (CO 2 / 02/ N 2 ) across blood-brain 

barrier 

Volumefric fraction of j (CO 2 , 02/ N 2 ) 
in dry alveolar gas (i ~ A), 
in dry expired gas (i -E), 
in dry inspired gas (i = 1). 

Blood oxygen capacity 

conversion factor 

Dissociation constant for 
carbonic acid 


K| Volume of brain (i = B), 

cerebrospinal fluid (t - CSF), 
tissue (i = T), a! .^.oii (i = L) 

MR;j Metabolic rate of carbon dioxide 

production (j -CO 2 )/ 
oxygen consumption (j — O 2 ) by 
brain (i = B), tissue (5 — T ) 

Pj j Partial pressure of | (CO 2 , 02/ N2) 

in alveoli (i - A), in blood at 
lung exit (i = a), in brain (i ~ B), 
in cerebrospinal fluid (i -"'CSF), in 
inspired aii (i I). 


Units 

1 . liters x/liter y 
where x - gas 
or 02 for Hb 02 
and y -blood, 
brain, CSF, or 
tissue . 

2. nanomol es/li ter y 
For H 1 ion and 

y defined as above 

I iters/mi n per 
mm Hg 

dimensionless 


liters/liter blood 
atm/mm Hg 
nanomoles/liter 


liters 


iiters/min 


mm Hg 
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Symbol 


P H; 


Q 


Qi 


AQ 


AQ; 



t 



V; 


Description 


pH of blood at lung exit (i = a), 
of blood at brain exit (i = vB), 
of blood at tissue exit (i = vT), 
of cerebrospinal fluid (i =CSF) 

Cardiac output 

Cerebral blood flow (i ==B) 

Normal resting blood flow fi =N) 

Change in cardiac output 

Change in cerebral blood flow (i = B) 


Units 


liters/ min 
liters/min 

liters/ min 
litets/min 


Change in cardiac output due liters/min 

to CO 2 (j =C02), due to 

o 2 (i =o 2 ). 


Time constant associated with first order minutes 
differential equation for cardiac output 
response (5 = 1), for cerebral blood 
flow response (i = 2). 


Time 


minutes 


Blood transport delay 
from lung to brain (i = aB) 
from lung to tissue (i ~aT) 
from brain to lung (i ~ vB) 
from tissue to lung (i = vT) 
from lung to carotid body (i = ao) 
for specified path segment (j =(!)) 

Expiratory gas flow rate (i = E) 
Inspiratory gas flow rate (i = |) 


minutes 


liters/min 


One of the major modifications of the basic Grodins' model involves the addition 


of the exercise (workload) simulations. Table 2 gives a description of some of the terms 
used in the modification's subroutines. 
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Table 2. Symbols and Terms Associated With the Exercise Subroutines in the 
Modified Grodins 1 Respiratory Control System Simulation- 


Symbol 

Description 

Units 

CXT 

Term set equal to simulation time 

minutes 

DEADVT 

Dead space ventilation 

liters 

DU RAT 

Workload duration 

minutes 

FREQ 

Respiratory frequency 

minute ^ 

HRATE 

Heart rate 

minute ^ 

RMUN 

Linear function of O 2 consumption 
as used in subroutine RC12 

liters/min 

RMTB 

Function of steady -state O 2 
consumption at specific workload. 

liters/min 

RMTB(2) 

Term set equal to RMT (2) in subroutine 
RC12. 

liters/min 

RMTM 

Term set equal to RMT (2) In 
subroutine RC 12. 

1 iters/mi n 

RMT(l) 

Metabolic rate of CO 2 in tissues 

liters/min 

RMT (2) 

Metabolic rate of O 2 m tissues 

liters/min 

SS02W(WORK) 

Steady-state O 2 requirement for 
specific workload 

liters/min 

SSVENT 

Steady-state ventilation for the 
workload of inteiest 

liters/min 

SV NT 

Parameter used to signify start of 
linear change in ventilation to 
steady -state value 

liters/ min 

SV NT 2 

Difference in ventilation between Vi 
chemical and the steady-state VI foi 
a paiticular workload 

liieis/min 
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Symbol 

Description 

Uni ts 

TCT 

Inverse of time constant associated 

dimensi onless 


with exponential components in 
exercise subroutine calculations 


TIMEON 

Time when workload is initiated. It 

minutes 


is set equal to CXT after workload 
is read. 


TV NT 

Minute Volume 

liters/min 

WORK 

Workload 

watts 

WORK2 

Woikload (new) 

watts 

2.2 Terminology 

Associated with Mil horn's Models 


Generally speaking, many variables appearing in Grodins' basi 

c model also appear 

in Milhorn’s models with different symbol notations. Table 3 gives 

a listing of most of 

the variables that 

are incorporated into the discussions of Milhorn's 

models as presented 

in Section 4 and Appendix 6.2. 



Table 3. Symbols and Terms Used in Milhotn': 

s 


Respiratory System Models 


Symbol 

Description 

Uni ts 

A 

Effective diffusing area of pulmonary 
membrane 

2 

mm 

-A 

Normal effective diffusing area of pul- 

2 

mm 

o 

monary membrane 


A, B,C,D 

Empirical constants used in Lloyd and 

A, B,C (mm Hg) 


Cunningham's controller equation 

D (liters/min mm Hg) 

Ca co 2 

Alveolar CC> 2 concentration 

vol fraction 


Alveolar 0 2 concentration 

vol fraction 

C A.o s 

Carotid bodies' site 0 2 concentration 

vol fiaction 
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Symbol 

Description 

Units 

Cb co 2 

Brain tissue C 0 2 concentration 

vol fraction 

%' 

H + concentration at central receptor 
on blood side of blood-brain barrier 

nanomoles/liter 

C Bo 2 

Brain tissue C>2 concentration 

vol fraction 

C CSF^ 

H 1 concentration at central receptor 
on cerebrospinal fluid side of blood- 
brain barrier. 

nanomoles /1 iter 

C PS H t 

H + concentration at peripheral sensors 

nanomoles/liter 

%o 2 

Body tissue CO2 concentration 

vol fraction 

C To 2 

Body tissue 0 2 concentration 

vol fraction 

D 

Brain tissue CO2 diffusion coefficient 

liters/min/mm Hg 

f 

Blood flow through a single capillary 

liters/min 

f V D 

Minute dead space 

1 iters/mi n 

k; ; 5=1,... ,8 

Constants used in cerebral circulatory 
controller equation 

— 

k 2 

Constant determined from CO2-O2 
dissociation curve 

dimensionless 

k 5 

Ratio of normal arterial PO2 *° normal 
alveolar P02 

dimensionless 

n 

Constant used in a controller equation 

dimensionless 

P 

a C 0 2 

CO2 arterial tension 

mm Hg 

p 

a o 2 

C>2 arterial tension 

mm Hg 

Pd 

“CO 2 

Deep brain tissue CO2 tension 

mm Hg 

(> csfco 2 

Cerebrospinal fluid CC>2 tension 

mm Hg 
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Symbol 

Description 

Units 

p h 2 o 

Water vapar pressure 

mm Hg 

P|C0 2 

Inspired carbon dioxide partial pressure 

mm Hg 

p l° 2 

Inspired oxygen partial pressure 

mm Hg 

6s 

Cerebral biood flow 

iiters/min 

Qc 

Total blood flow through pulmonary 
capi Maries 

liters/min 

Q,/~ 

v -o 

Normal blood flow through pulmonary 
capi 1 laries 

Iiters/min 

^CS 

Average blood flow at effective sensor 
si te 

liters/min 

Q S 

Pulmonary shunt flow 

liters/min 

R 

Respiratory quotient 

dimensionless 

s 

CO 2 dissociation curve's slope 

mm Hg ^ 

e 

V 

Minute pulmonary ventilation 

iiters/min 


Alveolar ventilation 

liters/ min 

v co 

Rate of change in volume due to C0 2 

liters/min 

2 

production 


V D 

Dead space 

liters 

V L 

Lung volume 

1 i tei s 

V L 

n 

Normal lung volume 

liters 

o 

V Oo 

Rate of change in volurhe due to 0 2 

liters/min 

2 

consumption 


V T 

Tidal volume 

liters 

X 

Effective depth of sensor below 
medulla's suiface 

mm 
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3. EVALUATION OF G RODINS' RESPIRATORY . 

CONTROL MODEL MODIFICATIONS 

3. 1 Statement of Objectives 

One' phase of the research effort was the evaluation of Grodins' respiratory control model 
as it simulated- the effects of altered environments; in particular, the simulation of various 
exercise levels in an altered environment. Part of this research is described in a Study Report 
(See Appendix 6. 1). Additional investigation involving the interactions of exercise and 
respiratory control is included in this section. 

The objectives of this component of the research was to evaluate Grodins' program and to 
recommend modifications to increase the program's fidelity. This included the establishment of 
an evaluation procedure, comparison of program responses to published data, identify problem 
areas, and indicate modifications to rectify or alter known simulation shortcomings. 

3.2 Description of Grodins' Respiratory Control Model and Modifications 

The approach used to present the material in this section will be in the form of a 
summary of the Study Report (Appendix 6. 1). The portions of the Study Report that were 
expanded after it was written, i.e. further investigations involving the exercise subroutines, 
will be presented in more detail here. The system can be envisioned as a feedback control 
system comprised of a "plant" (the controlled system) and the "regulatory component" 
(controlling system). 

3.2.1 Controlled System 

The controlled system as represented in Appendix6. 1 Is partitioned into three compartments 
corresponding to the lungs, brain, and tissue with a fluid interconnecting path representing 
the blood. A set of differential-difference equations describing a gas transport and exchange 
system between the lungs, blood, brain, and tissue compartments constitutes the framework 
for the model. Equilibria formulations involving acid-base buffering and metabolism for the 
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alveolar-arterial , venous blood-brain, and venous blood-tissue inteifaces are qualitatively 
and quantitatively described. As indicated these formulations are the most complex of any 
of the model segments. 

Reference should be made to Appendix 6.1 for a condensed listing of the assumptions and 
features of the model's compartments. It is important to note these assumptions since an 
evaluation of the system's performance would be irrelevant without them. 

An important aspect of the simulation is related to the blood transport time delays. 

As with any physically dynamic system involving flow dynamics or transfer of information, 
there are time delays in signal transmission. The time delay, 'C, is associated with the blood 
flow from the lung to the brain compartment, from the lung to the tissue compartment, from 
the tissue to the lung compartment, from the brain to the lung compartment, and from the 
lung compartment to the carotid body receptors. These are very necessary terms and their 
incorporation in the system equations enhances the confidence placed upon the dynamics 
of the simulation. To further define these time delays, they are not considered to be 
predetermined constants but are continuously recalculated. They are time dependent in 
that they are functions of cardiac output, brain blood flow, and vascular segment volumes. 

This reasoning is very sound since the blood flow dynamics obviously affect relative site 
gas concentrations. It should also be mentioned that the inclusion of these variable time 
delays necessarily increases the time for the calculation routines. Therefore, the sophistication 
that they provide might be offset by the cost of the simulation process. 

3.2.2 Controlling System 

To provide the closed loop (feedback contiol) phenomenon a controlling system is 
utilized. Ignoring the exercise influence, thiee of the controlled system's outputs aie 
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monitored as inputs for the controlling system. A controller equation describing the inspired 
ventilation as a weighted function of the H + concentration in the cerebrospinal fluid 


concentration 
and the oxygen contribution at the carotid 


^CcSF(H + )^' ^e concentration in the arterial system ^(H+^aoj)' fhe H 
in the venous blood of the brain t-) (t ~ T 0 $)J r 

chemoreceptor sites (F A (q 2 }) provides this feedback mechanism. Using the terminology as 
defined in Appendix 6.1 the controller equation is written as 


V, - 1.1380 [C(16)C a(H+) (t-r< aB ) h (1.0 - C(16) ) C CSF(H -t } )] 

+1.1 540^+} (t -^q) -f TERM - VI (N) (3-1) 

where 

TERM = 23.6(10)" 9 [104 - (B-47)F a( q 2) ( t - ^ a0 )] 4 * 9 

for (B - 47)F A(Q2) (t < 104 (3-2) 

TERM = 0 for (B - 47)F A(c , 2) (t -Y QO ) > 104. . (3-3) 

The study has indicated that the controller equation is the most flexible component of 
the entire simulation and also the component which provides the most uncertainties. This fact 
is even more evident when variations in exercise levels are being simulated. Equation 3-1 
takes on a different form when exercise levels are being simulated. Several other forms of 
controller equations are presented in both Apprendices 6. 1 and 6.2. The idea that the controller 
equation needs continued evaluation is amplified in the conclusions and recommendations. 

3.2.3 Exercise and Related Subroutines 

in Appendix6.1 several simulation runs were made in order to establish the significance 
of the exercise subroutine-modification that Weissrrlan (28) had added to Grodins’ program. 

The ventilation responses foi varying exercise levels at an alteied enviionment (Skylab 
conditions) did not appear to be completely satisfactory. Fiist of all, the neural component 
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response in Hie off -transient was not evident. Also, after a change in the exercise level 
was initiated the transient response to a steady-state value of ventilation seemed to be 
longer than would be desirable. To better understand the exercise subroutines a more 
detailed investigation of the computer progiam was performed (27). The following material 
is a portion of that investigation. The existing computer program (Appendix 6.3) has a limited 
number of comment statements. Fortunately, the entire program is partitioned into subroutines. 
The exercise modification is contained in subroutine RC12. The block diagram in Figure 3 
shows subroutine RC12's relationship to the other subroutines in the program. RC12 is 
embedded in a loop that establishes and solves the controlled and controlling system equations. 
Therefore, the -control I er subroutine, RC17, is also a part of the loop. The controller 
subroutine is called before the system equations are calculated in RC11 since inspiratory 
ventilation rate (Vi) and expiratory ventilation rate (VE) are both needed in the calculation 
of the system equations. RC17 will be considered later. 

A flow chart shown in Figure 4 was formed in order to get an insight into the functioning 
of the exercise subroutine. In this subroutine the tissue metabolic rates of oxygen (RMT{2)) 
and carbon dioxide (RMT(1)) are calculated with heart rate (HRATE), respiratory 

frequency (FREQ), minute volume (TVNT), and dead space (DEADVT). The above mentioned 
metabolic rates are the variables that interface this subroutine with the remainder of the 
program. Both metabolic rates are used in the system equations. 

The subroutine reads in a workload (WORK2) and duration (DURAT) from the input data 
cards. See Appendices 6. 1 and 6.3 for a description of input data cards. A decision is then 
made concerning the workload which depends upon the workload's absolute magnitude and 
its magnitude relative to the previous workload (WORK), This decision is needed since theie 
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FIGURE 4. Flow chart for RC12 
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FIGURE 4(c). Continuation of Flow chart for RC12. 
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FIGURE 4(d). Continuation of flow chart for RC12. 
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FIGURE 4(f). Continuation of flow chart for RC12 









30 



Fig 4 (U) 


FIGURE 4(g). Continuation of flow chart for RC12. 
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FIGURE 4(h). Continuation of flow chart for RC12 
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are two different equations that calculate the metabolic rate of oxygen. One is for WORK2 
less than WORK, and one is for WORK2 greater than equal to WORK. To make sure 
that the program knows where it is at all times, a flag is set depending on the relative 
magnitude of the new workload. If WORK2 is less than WORK, MARKER is set equal to 
zero (MARKER '• 0). If it is greater than or equal to WORK, MARKER is set equal to one 
(MARKER - 1). This flag is necessary because the program is continually calling RC12. 

The flag makes it possible for the program to return to the same set of equations until a 
new workload is read in. This will occur whenever the simulation time, C(3 5), is greater 
than TIMEOF. TIMEOF is set equal to the duratioh time plus the simulation time, CXT, when 
the workload is read in: 

TIMEOF = CXT + DU RAT (3-4) 

Physiological studies show that oxygen uptake is a good indicator for the amount of work 
that is being performed (31). For this reason many other physiological parameters are 
always related to oxygen uptake. As mentioned before, the modification under discussion 
for this program is based mainly on two expressions for oxygen uptake; one for when the 
workloads are increasing, and one for when the workloads are decieasing. These two 
expressions are discussed in the following material. 

(A). Increasing Workloads 

'For this case subroutine RC12 utilizes the expression 

RMT(2) = SS02W(WORK) - (SS02W(WORK) - RMTB2)* 

EXP(-TCT (CXT - TIMEON)) (3 - 5) 

where: 

SS02W(WORK) is the tissue steady-state oxygen consumption for a particular woikload. 
Oxygen consumption varies linearly with workload as shown in Figure 5. However, the 
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relationship between oxygen consumption and workload in this program is a piecewise 
linear relationship, as illustrated by Figure 6. Figures 5 and 6 illustrate a maximum value 
(3.5 iiters/min) for the oxygen consumption. At this point there is no additional increase 
In oxygen consumption for a further increase in workload. This value is known as the maximum 
oxygen uptake (MAX O 2 / It is a measure of one's physical fitness. In the average male 
the MAX C>2 is about 3.5 liters per minute at a workload of 250 watts. MAX C >2 is higher 
in the more physically fit individual. The dotted line in Figure 6 is a suggested change that 
would eliminate the discontinuity at 250 watts. 

RMTB2 is the oxygen consumption at the time the workload is initiated. If the workload 
is started from rest, this value is 0.215 liters per minute. 

TCT is the inverse of the time constant. It determines the rise time of the exponential 
components. TCT is related to the workload by the following expressions: 

4.6 for work ^50 watts 

TCT = 

2 . 3/2* WORK/200 for work > 50 watts (3 - 6) 

The term TCT as a function of workload is plotted in Figure 7. From this plot it is seen 
that, as the workload is increased, TCT becomes smaller. This means that as the workload 
increases, more time is required to reach steady-state . 

CXT is equal to the simulation time C(35). 

TIMEON is the time when the workload is started. After the workload is read in, 
the subroutine sets TIMEON - CXT. Thus CXT - TIMEON is elapsed time of the workload 
at any time C(35). 

To gain a better insight into the functional relationship of RMT(2) as workloads are 
varied. Figures 8 and 9 illustrate the times required for the steady-state oxygen uptakes 
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to be achieved when exercise transition is from rest or from other exercise levels. In both 
figures it is readily seen that the transient response is longer for greater workloads, 
the transition being either from rest or another workload. 

(B). Decreasing Workloads 

For decreasing workloads RC12 utilizes an oxygen consumption expression that 
is a decaying exponential. It has a time constant that is twice as long as the time constant 
used in the expression for increasing workloads. The expression is written as 

RMT (2) = RMTB - (RMTB-RMTM) * EXP(-t CT*(CXT-TlMEON) J 0.5) (3 - 6) 

where: 

RMTB is the steady-stqte oxygen consumption at the new, lower workload. 

CXI, TCT, and T1MEON are defined as with increasing workloads. However, it is 
important to notice that TCT is determined not by the new workload but by the previous 
workload. This is physiologically justified because the swiftness by which the body returns 
to the resting state is dependent upon the prior exercise level. The greater the prior exercise 
level the longer if takes to replace the oxygen and high energy stores of the body. 

Figure 10 is a plot of RMT (2) for transitions from various exercise levels to rest, for 
five different workloads. 

No changes are made for the metabolic production of CC >2 during increasing and decreasing 
workloads. For both increasing and decreasing workloads the metabolic production of 
carbon dioxide, RMT(l), is calculated by the same equations: 


RMT(l) -- (TV NT i 40 . 77 )' RMT (2)/88 . 5 
for TV NT less than or equal to 37. 


(3-7) 





For TVNT greater than 37 

RMT (1 ) = 0.88RMT(2) 

where TVNT is the minute volume. 
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(3-8) 


(C). System Controller Equation 

Up to this point, only the alteration of the controlled system's equations during 
exercise has been discussed. It is still necessary to determine how the controller equation 
subroutine, RC17, is modified to take exercise into account. 

As before, a flow chart can be obtained for RC17 as shown in Figure 11. From the flow 
chart it is still difficult to determine the interaction of the steps and how the various 
individual elements of the routine interact. A print routine which prints out the important 
variables is added to this subroutine. See Appendix 6.3 (RC17) . When the program is run 
at different workloads, knowledge of the subroutine functioning is obtained by utilizing these 
printouts. 

Verification that the modification of the controller equation is based on the neurohumoral 
theory as discussed in Section 1.3 is possible. That is, it consists of a fast neural component 
and a slower humoral component. Referring to Appendix 6.3, which is the subroutine 
listing, the first calculation made is that of the ventilation as given by Gredins' controller 
equation, (VI). This value is then subtracted from the steady-state ventilation for the workload 
of interest, (SSVENT(SS02W(WORK))), This steady-state ventilation is calculated using 
the relationship given by the plot in Figure 12. This plot effectively illustrates that steady- 
state ventilation is linearly related to steady-state oxygen consumption. However, oxygen 
consumption is dependent upon the workload as pieviously discussed. The difference. 

If less than or equal to 15 liters per minute, between' this steady-state value and the value 
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obtained by Grodins 1 controller equation is the neural component. If the difference is 
greater than 15 liters per minute, the neural component is set equal to 15. The difference 
is given the variable name SVNT2. This neural component is added to the value of ventilation 
that is calculated by Grodins' controller equation and a new value for VI is obtained. 

Thus VI is defined as 

VI VI i SVNT2 foi 0 * SVNT2 15 (3-8) 

or 

VI ~VI + 15 for SVNT2 >15 (3-9) 

At this point, another variable is calculated, SVNT. SVNT is the difference between 
another steady-state value for ventilation defined by the term SSV ENT(RMLIN) and the 
value of ventilation which Is calculated in Equation 3-8 or 3-9. The variable SVNT is 
formulated as 

SVNT = SSVENT(RMLIN) - VI . (3 - 10) 

The argument RMLIN used to determine the steady-state ventilation term Is calculated in 
RC12. RMLIN is a slowed down, [Inearlized version of RMT(2). It is given by 

RMLIN =SS02W(WORK) - (SS02W(WORK) - RMTB2)(1 . VTIME) (3-11) 

where SS02W(WORK) and RMTB2 are as previously discussed. VTIME is the function that 
slows down RMLIN in relationship to RMT(2). It is given by 

VTIME =TCT(CXT -TIMEON)/9.2 for VTIME <1 (3-12) 

The value 9.2 makes RMLIN approach its steady-state value at a rate four times slower than 
RMT(2) when the workload is 100 wafts. For workloads less than 100 watts, the response of 
RMT(2) is less than four times as fast. For woikloads greater than 100 watts, the response 
of RMT (2 1 is greater than four times as fast Figures 13 and 14 show the relationship between 
RMLIN and RMT (2) for two different workloads (25 watts and 75 watts). Figure 13 represents 
a hansition from rest to 25 watts while Figuie 14 represents a transition from 25 wafts to 75 watts. 



OXYGEN UPTAKE (L/MIN) 



FIGURE 13 . The relationship between oxygen uptake and RMLIPjS 
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FIGURE 14. The relationship between oxygen uptake and RMLIN 
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When VTIME reaches the value of 1, RMLiN has reached steady-state and Is given by: 

RMLIN = SS02W(WORK). (3-1?) 

Returning to the discussion of the controller equation, if SVNT is greater than 0.5, 
ventilation will linearly approach its steady-state value. This is done by adding 75% of 
SVNT to the value of VI calculated above. The relationship used is: 

VI = VI + 0.75SVNT ' (3-14) 

After inspiratory ventilation has been calculated, expiratory ventilation is calculated, 
again using a relationship given by Grodins (2). 

Figure 15 illustrates plots of the various variables discussed above for a workload of 
75 watts. From these plots it is easily shown that the sum of VO and SVNT2 equals SSVENT 
(SS02W(WORK)) for SSVENT(RMLIN<0.5 + VI. VI is equal to the sum of VO and the 

> 

neural component. When SSVENT (RMLIN) becomes greater than the sum of 0.5 + VI, VI 
starts a linear increase to its steady-state value for that particular workload by using Equation 
3-14. 

The controller equation modification has no physiological significance. This is justified 
by the detailed investigation of the appropriate subroutines. As indicated by Strippoli with 
further justification presented in Section 3.4 it does permit a reasonable simulation for 
specific workloads (27). The fast on- and off-transient responses can be associated with 
neural components. However, the functional expressions which simulate this component do 
not actually describe the physiology involved. This particular physiology has not been well 
documented. 

3.3 Model Parameter Variations 

One of the goals of the investigation was to evaluate Grodins' modified program with 
regard to sensitivities for reasonable variations in parameter values. Grodins' modified 




TIME(MIN) 

FIGURE 15. Relationship between important variables used in 
controller equation at 75 watts 
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program provides for various parameter variations by the use of input data cards. There <?re 
many standard conditions for inspired volumetric fractions of gases from which to base a 
sensitivity analysis. The decision was made to use the following data from which all 
simulations of variations in parameters are compared 

B - 760 mm Hg 
FI(C02) = 0.1000 
Fl(02) - 0.1100 

FI(N2) = 0.7900 

See Appendix 6.1 for the type of responses obtained and the table of input data cards 
used for the base run. 

It is believed that these base conditions were adequate for the sensitivity analyses 
performed. However, one should use extreme caution in projecting, in a direct fashion, 
the sensitivity of any parameter to other base conditions. Only one parameter or group of 
common parameters was altered at one time. This is really the first step in what should be 
a continuing effort in analyzing combinations of parameter variaiions and weightings. 

The format followed in the Study Report (Appendix 6.1) is one which hopefully conveys 
the model's capabilities and limitations in a concise and informative manner. There is 
a tremendous amount of data which is not discussed or presented in graphical form. Only the 
variables which present significant changes are addressed. This is not meant to deemphasize 
the other variables' importance. General Electric Company's Technical Information Release 
741-MED-3008 and 74I-MED-3009 illustrate the output data format and the possible relationships 
between variables that can be generated (32,33). 

With regard to inspired gaseous concentrations only the on-transient and related steady-state 
conditions are analyzed. The off-transient as related to the transition from an exercise to 
the resting state is presented later in Section 3.3,2. 
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3.3.1 Resting State and Sea Level Conditions 

The computer program is a modified version of Grodins 1 program similar to that used 
by Weissman (28). Although it has a subroutine which involves a work-exercise phenomenon, 
the case of the resting subject is considered first. A table in Appendix 6. 1 presents the 
input data format giving card number, symbol, normal initial value, and a brief description 
of each parameter. 

Since a detailed discussion and evaluation of the parameter sensitivity analyses are 
presented in Appendix 6.1 only a few comments are included here. 

The controlled system parameters that were varied include 

(1) The time constants R1 and R2 of the first order differential equations 
describing cardiac output and brain blood flow, 

(2) the standard bicarbonate content BHC03 BRAIN and BHC03 TISSUE 
of the brain and tissue compartments, and 

(3) the diffusion coefficients D(C02) and D(02) of CO 2 and O 2 across the 
blood-cerebrospinal fluid barrier. 

Briefly, the simulation proved to be insensitive to variations in R1 and R2 after an initial 
transient disturbance of about 1 minute duration. The simulation was found to be much more 
sensitive to changes in the bicarbonate content of the tissues than of the brain which is thought 
to be related to relative volume sizes of the two compartments in addition to the buffering 
and metabolism relationships. Results of the simuldtions with varying diffusion coefficients 
supported the known physiology (blood-cerebrospinal fluid barrier highly permeable to CO 2 ) 
involved in that the system was sensitive to D(C02) variations and insensitive to 0(02) 
vaiiations (34). 
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The controlling system parameter weightings that were varied for the resting subject 
include 

(1) the weighting of the H + concentration in the CSF compartment, CNT SENS 
COF, 

(2) the weighting of the H + concentration at the carotid bodies' site, CRTD 
BDY SCF, 

(3) the weighting of the relative influence of the H + concentration in the CSF 
and brain compartments, CENT SENS PT, and 

(4) combinations of (1) and (2). 

No significant variations were observed in any compartmental FT 1 " concentration, 

P 

a 

directly dependent upon CNT SENS COF although not responding rapidly in the initial 
transient. Considering all of the variables monitored there was a greater degree of sensitivity 
toward variations in CRTD BDY SCF than CNT SENS COF. Thus, there is an indication 
that the system responds more rapidly to chemical changes at the carotid bodies' site than 
in the CSF compartment. This feature was further emphasized when combinations of 
CNT SENS COF and CRTD BDY SCF variations were used. Simulation runs illustrated that 
the brain compartment H + concentration detector mechanism responds faster than the CSF H + 
concentration defector system. All of the variables illustrated a faster response as CENT 
SENS PT was increased. 

3.3.2 Exercise States at an Altered Environment 

The usefulness of such a model must be evident for conditions which approximates those 
of actual experiments. Thus, this model must perform quality simulations under environmental 


(O 2 ), ° r ^ a (C02) ^ or varia, '' ons ' n CNT SENS COF. However, VI was shown to be 
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conditions such as those of Skylab. From the simulations run for sea level conditions, 
barometric pressure -760 mm Hg, it was concluded that the most sensitive component of 
the system was the controlling equation. Thus, the simulation efforts related to an altered 
environment were concentrated on the controlling equation and its variable weightings both 
for the resting state and various exercise levels. 


Initial variable values were obtained by allowing the simulation to reach steady-state 
conditions when subjected to the altered environment. This first simulation effort involved 
a change of only the barometric pressure and inspired gaseous concentrations. Thus, 
these specific input data cards were altered as shown in Table 4 with the CC >2 concentration 
corresponding to 5 torr. 

Table 4. Input Data Cards Specifying the Altered Environment 


Card No. 

Symbol 

Revised Initial Values 

C(30) 

B 

260.0000 

C(31) 

FI(C02) 

.0192 

C(32) 

FI (02) 

.7000 

C(33) 

Fi(N2) 

.2808 


The use of the steady-state values obtained from this run compensated for the preconditioning 
period. See Table 2 of Appendix 6.1. Exercise levels of 5-minute intervals were utilized 
as shown in Table 5 for the following set of time intervals. 

Table 5. input Data for Workload Schedule 


Time Interval 


Workload * 


0 < t < 5 min 0 watts 

5.min <t< 10 min 50 watts 

10 min <t< 15 mtn 100 watts 

15 min <t< 20 min 150 watts 

20 min 30 min 0 watts 


*See Amend! x 6.3 for Workload Input Data format. 
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Since the simulations discussed in Section 3.3.1 indicated that the most influential 
component of the system was the controller equation simulations of exercise at Skylab 
conditions concentrated on further evaluations of the controller equation. 

The controlling system parameter weightings that were varied during the transition 
from one exercise level to another include 

(1) the weighting of the relative influence of the H + concentration in the CSF 
and brain compartments, CENT SENS PT, 

(2) the weighting of the H + concentration in the CSF compartment, CNT SENS 
COF, and 

(3) the weighting of the H + concentration at the carotid bodies' site, CRTD BDY SCF. 

All exercise levels produce a slightly lower steady-state response for V| when CENT 

SENS PT is increased while the decrease in exercise level is not significantly altered by 
variations in CENT SENS PT. The greater the weighting of the H + concentration sensing 
in the venous blood of the brain, the higher the relative arterial CC>2 tension is for any 
instant of time. As expected the level of the H + concentration in the CSF compartment is 
greater when CENT SENS PT is increased. 

Since Ccsp(H + ) does not deviate from the base run for variations in CNT SENS COF 
other variables are not altered significantly either. As observed in Section 3.3. 1 the 
simulation is more sensitive to alterations in CRTD BDY SCF than CNT SENS COF. This fact 
is also evident under exercise conditions. 

Appendix 6.1 also includes plots of P a (Q 2 ) versus Q for ^e complete exercise cycle 
(rest-exercise levels-rest). One interpretation of these plots is that a depletion and lecovery 
of arterial oxygen is evident for the exercise transitions that do not initiate from a resting state. 
Other conclusive statements should not be made without further investigation. 



55 


3.4 Evaluation of Exercise Feature 

For further supportive evidence of the role that exercise plays In the control of respiration, 
other exercise levels and level transitions were simulated (27). In particular, the justification 
of both neural and humoral components was considered with an expansion of the investigation 
of the exercise subroutines discussed in Section 3.2.3. The values for the input data cards 
are shown in the following listing. All other input data cards hove values indicated in 
Appendix 6. 1 . 


Table 6. Input Data Cards for Five Exercise Simulation Runs 


Input 

Data 

Cards 

Run 

Run ^2 

Run ^3 

Run H 

Pun ^5 

C(l) 

0.1783 

J 

r 

0. 1674 

0. 1746 

j 

\ 

C(2) 

0.5336 



0,5419 

0.523 1 



C(3) 

0.2881 


j 

0.2907 

0.3123 



C(4) 

0.6413 



0.6345 

0.6359 



C(5) 

0.0012 



0.0011 

0,0011 



C(6) 

0.0011 



0.0011 

0.0011 



C(7) 

0.6153 

Some 

0.6164 

0.6390 

Same 

C(8) 

0.0015 

as 

0.0011 

0.0007 

a? 

C(9) 

0.0012 

Run ^1 

0.0012 

0.0011 

Run #\ 

COO ) 

6 ,0000 



7.3800 

11.5717 

1 


C(ll) 

0.7496 



0.6839 

0.7260 

j 


C(12) 

48. 1202 



47.9542 

47.3567 



C<13) 

36.6316 

i 


36.0486 

35.3290 



C(I4) 

70.6804 

1 


68.5268 

66.2437 



PMT(J) 

0. 1820 



0. 1820 

1 . 1699 



RMT(2) 

0.2150 

V 

t 

0.2150 

1.3294 

\ 

/ 

WORK2 

0.0/0. 2 * 

0.0/0. 2 15 

25/0.2 ' 

100/2.5 * 

0.0/0. 2 * 

WORK2 

50/1.8 * 

25/1.8 * 

70/2.0 * 

0.0/30 * 

50/2,5 - 

WORK2 

100/2.0 * 

75/2.0 * 

100/2.0 * 



0. 0/2.0 * 


*Watts/MIn 

Gallagher (22) showed that at workloads of 50, 100, and 150 watts no neural component 


was simulated for the transitions to 100 and 150 watts when the tiansition was made duiing 
steady-state. The absence of the neuroi response was also shown by some of the above 
simulation runs. Figure 16 is a plot of VI for Run ^1. 
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For a simulation run of 0, 25, and 75 watts (Run ^2) neural components were observed 
for both exercise transitions. Figure 17 illustrates the VI response for this run. It is also 
observed that the neural response accounted for the entire response in the 0 to 25 watt 
transition. This type of response was also justified by D'Angelo et al. (23) when they reported 
that ventilation during low levels of exercise was piedominately associated with a neural 
response . 

Run # 3, having exercise transitions of 25 to 7D watts and 70 to 100 watts, simulated 
neural components at both transitions as shown by Figure 18. Although no quantitative 
experimental data is available for comparison it appears that the neural component is sensitive 
to the level and increment of exercise levels. Also, the problem discussed later in this 
paragraph, that of system initialization, prevailed in the initial transient response; however, 
the initial deviation in the system’s variables are not as severe for 25 watts. Run H, with 
exercise levels of 100 and 0 watts, was used to determine the type of response that would be 
obtained in the off-transient of exercise. The off-transient response obtained was not 
justifiable. Further investigation of the exercise subroutines revealed that at the first 
calculation of the controller equation, which occurs before RC12 is implemented, the 
initial work level is considered to be 0 watts instead of the desired 100 watts. Thus, the 
transient responses are invalid. Further modification of the program is required if simulations 
are to be initialized with other than 0 workload (or small workloads). 





V i (L/MI N] 



TIME(MIN) 

FIGURE 18. Simulation run number three for 25-70-10 0- 0 w a 1 1 s 
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4. MILHORN'S RESPIRATORY CONTROL MODELS 

4.1 Statement of Objectives 

! 

The other phase of the research effort was the study and analysis of Milhorn's respiratory 
control models. The major thrust of this effort was in the summa ization of the theoiy and 
assumptions of existing models which are associated with Milhorn's research (4, 5, 7, 8 , 24 , 35). 
One aspect of the objectives was to analyze published piogiam lesulls and describe potential 
applications of Milhorn's work to improve Grodins 1 piogram. These potential modifications 
could then be implemented in Grodins 1 program and evaluated at a future date if the con- 
tinuation of the modelling effort should include such implementation. 

4.2 Models' Variations and Significant Features 

The detailed evaluation of Milhorn's models leading to the fulfillment of the above 
objectives is presented in Appendix 6.2. Most of Milhorn's models, although developed 
under different experimental conditions and for varying motivations, contain features that 
are similar to Grodins 1 work. Therefore, to complete an investigation of two of the major 
groups of respiratory control models it seemed appropriate to discuss the salient features 
of Milhorn's models and suggest special modifications which might enhance Giodins' model. 


With Appendix 6.2 (Study Report) available it is not necessary to go into any depth 
in the summary of the Study Repoit. Only a brief statement concerning the vaiious models 

that were presented in the Study Report is contained in this section. 

« 

One of Milhorn's first models was a compai tmental lespuatory model with the basic 


assumption being that the control of alveolai ventilation was due to CO2 and O2 levels (7). 
The effects of the H* stimulus was related to 8(302 ^7 a threat relationship. The model 
contained a controlled system and a controlling system; thus, having an oveiatl similarity 
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to that of Grodins' model. Upon the establishment of variable arid functional interrelationships 
a set of eight variables and eight differential equations was used to describe the system. 

The ventilatory controller equation provided the feedback mechanism. An evaluation of 
some of the simulation responses is contained in Appendix 6.2. 

Another model which illustrated the static, dynamic, mechanical, and cyclic 
phenomena emphasized a restructuring of the ventilatory controller equation (8). Special 
weightings of the H" 1 contribution from the peripheral sensors, a central receptor in the 
brain, and a central receptor in the cerebrospinal fluid were given for the controller 
equation. No exercise phenomena were presented. The interrelationships between the 
cardiovascular and respiratory systems were experimentally investigated with the goal 


being to expand the range of Lloyd and Cunningham's controller equation range for 

The concept of venous admixture and a detailed description of the capillary shunts 


form the basis for another model (24). 


Features of this model would be of value to the 


present research if a modification of Giodins' model included a description of the gas 
exchange mechanism associated with the lung compartment. The most important advantage 
in the use of this model would be the component of the simulation associated with the 
admixture phenomenon as the fractional gaseous concentrations of the environment were 
altered. 


For a complete understanding of the transient ventilatory response, the receptor sites, 
number of sites, and the general nature of the responses at these sites should be better 
defined. Consequently, the system presented by Milhorn and Brown (5) was formulated 


in such a manner as to moie fully utilize the information horn the steady-state lesponses in 
oidei to impiove the tiansient simulations. Alveolai Oy and COy tensions aie descnbud in 
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terms of inspired O 2 , O 2 consumption, alveolar ventilation, respiratory quotient, water 
vapor pressure, and barometric pressure. Since all of these terms are accessible in Grodins 1 
model these equations could be made compatible with such a model. Interdependence 
of these forenamed variables are discussed in detail in Appendix 6.2. A comparison of 
Gray's controller equation and Lloyd and Cunningham's controller equation was made 
using the model framework as described above. Three-dimensional plots with inspired O 2 
and CC >2 tensions and either alveolar ventilation, tidal volume, or respirafoiy frequency 
were presented. Also three-dimensional plots with alveolar O 2 and CO 2 tensions and either 
alveolar ventilation, tidal volume, or respiratory frequency were given. With an expansion 
of this system to include exercise it appears to contain features which would be appropriate 
for suggested modifications. to Grodins 1 program. Further discussion of this modification is 
presented in Section 5. 

Stimuli of 3, 5, 6, and 7% CO 2 ~ air mixtures were used in a model developed by 
Reynolds et al . (35). The basic purpose of the experiment was to correlate the strength 
of the stimuli and the magnitude of the responses. Minute ventilation, tidal volume, 
respiratory frequency, alveoiar oxygen tension, and alveolar caibon dioxide tensions 
were recorded in order to better define the transient responses of the variables. Summary 
discussions concerning the transient half-times in addition to the responses of the 
forenamed variables for the various stimuli are presented in Appendix 6.2. This type of 
experimentation will help to chart subject responses for rapid changes in fhe gaseous mixtures 
of fhe air. 

The H + concentration is a variable that is well established as a contributor to the control 
of ventilation. The precise H" 1 concentration detectoi sites are not established. It is 
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conceivable that- two detector sites, one corresponding to the C5F compartment and one 
corresponding to the venous blood of the brain, might be appropriate for respiratory control 
monitoring of the H 1 concentration. A model has been formulated which attempts to 
establish a satisfactoiy transient response in ventilation which would apply if the CO 2 

1 

detection site (or H* concentration ) was assumed to be on either side of the control leceptors 
(4). With this model ventilation is assumed to be controlled by the peripheral sensor H 1 
concentration and the central sensor H h concentration. Pulmonary ventilation is the sum 
of two signals. These two signals have a weighted functional dependence upon the peripheral 
and central sensors. This special weighting allows for the establishment of a single sensor 
depth. Utilizing a CO 2 <ension gradient this depth was found to correspond to a position 
which is 77% of the distance between cerebrospinal fluid f*C02 an< ^ deep b ra m tissue Pq^- 
Three different controller equations were evaluated using this single sensor site depth. 

This section has contained a summary of a review of control system models developed 
by Milhorn and his associates. Chronologically speaking, the models encompass those 
developments associated with the steady-state models published in the mid-1960‘s to the 
more recent investigationsinvolving CFS perfusion studies. These recent developments 
also attempt to simulate transient responses as well as steady -srato responses. 
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5. CONCLUSIONS AND RECOMMENDATIONS - 
5.1 Grodins' Respiratory Control Model 

After an evaluation of the parameter sensitivity analyses, the simulation runs at 
an altered environment, and the exercise subroutines itappears that Grodins' Respiratory 
Control Model can be made adaptable to the research efforts of the Environmental Physiology 
Branch, Biomedical Research Division, Life Sciences' Directorate, NASA-JSC. Several 
alterations to the program which could improve the physiological soundness of the variable 
responses are included in this section. Here specific considerations are directed toward 
the exercise phenomenon while Section 5.2 contains suggested modifications for Grodins' 
program which evolve from Milhorn's work. 

With regard to the formulation of Grodins' basic piogram the parameter sensitivity 
analyses should be extended to include combinations of parameter variations especially 
in the off-transient environmental conditions. Considering the assumptions used in 
developing the model a tolerable range of environmental CO 2 levels should be established. 
Since levels of inspired CO 2 are physiologically correlated with ventilation rates, coma, 
anesthetized states, and death, the model should perform in a similar manner. If it doesn't 
perform as required CC >2 limits must be specified. 

At present, the model includes no general description of the control of cardiac output 
and regional blood flow. Certainly the legional blood flow feature will be important 
when the interfacing of the cardiovasculai , thermal, renal/endociine/body fluid, and 
tespiratory systems aie considered. This would pertain to both the testing and exeicismg 
subject . 

In the investigation of the intenelationships of paiameteis that compiise the exetcise 
subioutines of Giodins' modified piogtam it was obseived that the metabolic late of oxygen 
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consumption was the basis for the modification. The system controller equation which 
incorporates the exercise phenomenon is based upon the neurohumoral theory. Inspired 
ventilation (V|) is simulated by a fast neural component and a slow linear component when 
transition in exercise levels occurs. The slow linear component is related to oxygen con- 
sumption via a linear version of RMT(2) called RMLIN. See Section 3.2.3 and Appendix 6.3 
for the details of this modification. 

Based upon the investigations thus far, theie dre several suggested modifications for 
the existing program that would improve the simulation with regard to exercise. According 
to published experimental data, a transition in exercise levels produces a faster transient 
response for ventilation. With a modification of the expression for VT1ME (Subroutine RC12) 


from 




VTIME =TCT{CXT-TIMEON)/9.2 

(5-1) 

to 




VTIME -TCT (CXT-T!MEON)/4.6 

(5 - 2) 


the variable responses to exercise will be more acceptable. 

A punched card output is available for the first fourteen input data cards. If, in 
addition to these output data cards, the workload, time, metabolic rate of oxygen (RMT(2)), 
and the metabolic rate of carbon dioxide (RMT(l)) were included as punched output data 
it would permit the responses corresponding to one workload to be used as the initial 
conditions for a succeeding workload. To complete this change, the initialization of 
time and workload should be removed. Although not completely justified, these changes 
might also require a renaming of variables in the exercise portion of the controller equation 
since WORK is used for both resting and exeicisirtg conditions. 

The off-transient tesponse foi exeicise is not extiemely satisfying horn a physiological 
point of view. The output lesponse VI foi simulation Run ^5 (See the listing in Section 3.4. | 


is shown in Figuie 19. The initial off-transient contains an instantaneous deaiease in VI 
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followed by a linear decrease to the resting level. This instantaneous decrease which one 
associates with the neural response is of smaller magnitude than 'he corresponding instantaneous 
increase in the on-transient for 50 watts. Whether the implication that a neural component 
functions after exercise has ceased is physiologically sound or not is uncertain. 

it was observed that when an off -transient in exercise was initiated before steady-state 
conditions were reached for the previous exercise level there was no instantaneous decrease 
in the VI response. In fact, under these conditions, it was possible to generate an increase 
in the initial off -transient for VI. This particular response is caused by the fact that the 
routine for the off-transient is dependent upon the inspired ventilation reaching steady-state. 

In subroutine RC12, calculations are made for respiratory frequency (FREQ) and 
heart late (HRATE). These two functions are not utilized in any other portion of the program. 
They are given as functions of the C>2 metabolic rates of the tissue (RMBT(2)) and brain 
(RMB(02)) compartments. 

FREQ = 8.1+ 7.815 (RMT (2) + RMB(02)) (5 - 3) 

HRATE =43.8 (RMT,'?\ + RMB(02)) + 54.5 (5 - 4) 

These two expressions do not simulate in a faithful manner the respiratory frequency 
and heart rate for exercise transitions. For example, heart rate increases instantaneously 
in the initial transient of exercise and usually overshoots the steady-state value for any 
particular exercise level (31). But, since RMT(2j does not change in this rapid mannei , 
neither will HRATE. It is postulated that the immediate ipciease in heart rate corresponds 
to a neurological phenomenon. 

A fairly linear relationship has been established between steady-state heart rate 
and woth ioad. Staiting at appioximately 65-70 BPM foi the iesfing state the heait 
late mcteases to about 200 BPM foi maximal exetcise of a notmai individual. Data lelated 
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to these steady-state values in addition to data corresponding to an exercise dependent 
neurological component would form the basis for an improved HRATE expression. 

A similar postulation can be made for the respiratory frequency. In comparing the 
two formulations for FREQ and HRATE, FREQ has probably the more stable response of the 
two since the overshoot, if there is any, is not as pronounced as for HRATE. 

5.2 Mi I horn's Respiratory Control Models 

Several of Milhorn's models and their corresponding experimental justifications are 
presented in Appendix 6.2. The models that were summarized included those which 
concentrated on steady-state responses and also some of the more recent ones which simulate 
transient responses. The prime objective in the evaluation of these models was to establish 
some of their features as possible modifications for Grodins' model . As with Grodins' work 
several of Milhorn's models were involved in evaluation or establishment of controller 
equations. The suggested modifications involving the controller equation do not include 
any dependencies upon exercise. Thus, any controller equation alteration or substitution 
must be made compatible with the existing exercise controller functions of Grodins 1 model. 

The two most promising modifications of Grodins 1 program involve Equations 22 and 23 
of Appendix 6.2. They are repeated here for convenience. 


v - 2 (pa C o 2 " 37-24) (l + 13-6/^ - 25)) S 0 


r CS C 0 2 Pb C 0 2 r ^ r CSFco 2 ‘ B co 2 

Equation 5-5 allows for the additive and multiplicative effects of 0 2 and C0 2 instead 
of monitoring the compartmental H + concentrations. Since Equation 5-5 is only valid for 


+ ( P CSF,~r^ * P Br 


0 


exp -: 



Qcs s 
D/760 




(5-5) 


(5-6) 


P Ac0 2 ^ 37.24 mm Hg a new range is required so that simulations under the environmental 
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conditions of 5 torr CO 2 at 260 mm Hg can be run. Both variables involved here are 
accessible in Grodins 1 model. 

The modification involving the lone sensor site for the detection of CO 2 tension is 
more complex. This modification relates to Equation 5-6. The components of Grodins' 
controller equation which rely upon the H + concentration of the venous blood of the brain 
and the CSF subsystem would be replaced by a form of Equation 5-6. The essence of this 
modification is that it combines the sensing mechanisms of the brain and CSF compartments. 
Another technique involving the utilization of Equation 5-6 but still allowing the control 
of Vj to be H + concentration oriented would be to convert PcS^-q into a corresponding 
H + concentration. This particular H + concentration would be a weighted value corresponding 
to the effective sensor site. 

It is felt that these modifications would be the most beneficial to pursue in that they 


add missing detail and flexibility to the present modified version of Grodins' modei. 
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6. APPENDICES 

The Appendices Include two Study Reports and the program listing for the modified 
Grodins 1 Respiratory Control System. The Study Reports are essentially complete documents 
in that they each contain a table of contents, data with related discussions, conclusions, 
recommendations, and their own bibliographies. Both Study Reports and the program listing 
are included as Appendices in this research report to add continuity and completeness. 



Study Report - Respiratory Control System Simulation 
General Electric Company TIR-741-MED-3021 



f ^Tem R M 5 d?iI S - mmgrlzotion of fhe Mq i° r Features of Milhorn's Respiratory 


General Electric Company TIR-741 -MED-3030 



192 


6,3 _ Pr °9 ram Listing for the Modified Grodins 1 Respiratory Control System 

Minor modifications are indicated in the right hand margins of subroutines RC13 and 
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print routines as mentioned in the text. 
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INTRODUCTION 


I . 

A. Overview of Physiological System Modelling 

The importance of "biological system modelling within the 
area of biomedical research is becoming more evident. In 

i ' 

particular, the application of control theory which involves 
the mathematical description of element functionings and 
interrelationships between elements has made possible the 
• simulation of very complex physiological feedback control 
systems. The mathematical relationships range from algebraic 
formulations and ordinary linear differential equations to 

* 

nonlinear differential-difference equations and stochastic 

f !* ' * 

'''prbcesses , 

Prior to deciding upon a particular physiological system 
model it is instructive to utilize the expertise found in the 
literature and optimally meet the demands of the project with 
combinations or variations of existing research efforts. 

This is precisely the procedure established in this research 
involving the respiratory, system. 

B . Model Variations 

The system modelling effort that is discussed in the 
following report is one related to the respiratory control 
system. Several models have been developed for the respiratory 
control system emphasizing particular features which make 

J 

them adaptable for specific experimental and simulated 


environment a 1 condi Lions . 
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A few brief comments about some of the respiratory 
modelling efforts and their salient characteristics follow. 
Alyeolar ventilation responses were of prime importance for 
the output responses of a model developed by Bellville 
et al . (l). Frequency, amplitude, and mean level of inhaled 

C0 o were altered in a manner corresponding to sinusoidally 
varying CO g inhalation allowing a frequency response analysis 
of the respiratory 'system. 

I ' 

The phenomenon of pulmonary mechanics as related to 
various forms of inspired-expired flow systems has been 
researched by Fry .(2). Here variables related to physiological 
and pathological features of the air passageways are better 
defined. Elastic behavior of the lungs, airway compliances, 
and resistances were considered by Mead (3) in a model poss- 
essing characteristics similar to those of Fry’s model. 

Saidel et al. (it) developed a lumped parameter simulation 
of the pulmonary gas transport system with the ultimate 
goal of analyzing parameter changes. The model included 
five variable-volume compartments approximating the airway 
and alveolar regions and a const ant-volume compartment 
representing the capillary bed. Extensive efforts were taken 
mathematically describe the mass balance of the gas compart- 
ments and blood compartments. 

Breathing patterns related to COg oscillatory exchanges 

are presented in_a model developed by Yamamoto and Hori (5). 

- , 1 / » 

A linear controller equation is used to describe, indirectly, 

the dependence of tidal volume and rate of breathing upon 

I 



cerebral P 



(3) 
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Exercise and/or work load is a very important input for 
the proposed research. Analysis of maximum ventilation 
related to chronic airway obstruction and exercise has been 
conducted by Pierce et al. (6). This particular research was 
not intended to be a respiratory control system model study; 
however, the results of the endeavor are very applicable to 
modelling of airway passages and to the investigation of the 
respiratory-exerci se relationships . 

* Milhorn and Brown's steady-state human respiratory system 
(7) is very closely related to Gr.odins ' respiratory control 
system (8). In Milhorn ’'s model, a controlling system and a 
■controlled system form the basis from which mass balance 
equations are written describing the exchange mechanisms for 
'CO^ and O^. The controlling equation of Milhorn’s model is an 
important fundamental component which should be critically 
evaluated and compared to a similar component of Grodins 1 
model . 

Closely associated with Milhorn and Brown's model is a 
study of pulmonary capillary gas exchange by Milhorn and Pulley 
(9). This -work portrays a more detailed version of the capil- 
lary gas exchange by introducing anatomical parameters in the 
.mode'l. 

Assuming a second order representation for the mechanical 
p.ortion of the respiratory system Hilberman et al. (10) 
has utilised the Fourier series analysis and pliasor notation 
to determine the modulus of impedance, phase angle, compliance. 
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ry 

and resistance of the respiratory system. In doing so this 
research has developed a technique for adjusting parameters 
associated with respiratory mechanics. Although this research 
effort is not directly related to the present research the 
evaluation of parameters for specific environmental conditions 
might he useful "in extending the proposed model's capabilities. 
Since the pK of each compartment of the proposed model 

i 

is an important monitored variable the research by Hermansen 
et al. (ll) relating exercise and blood pH will prove to be 
a valuable piece of supportive research in determining a level 
of maximal exercise for the experimental subjects and in turn 
fo-r the simulation efforts. 

• A comparison of several respiratory system models' char- 
acteristics and capabilities are presented by Yamamoto and 
Raub (12). This review emphasizes that all models possess 
unique features but are categorized into three groups: (a) 
Models related to those of Grodins , (b) Models related to 

Morgan's efforts, and (c) Models related to Milhorn's research. 
1 \ 

A1e>o this review cites several other models giving their com- 
parisions to the above three groups. 

'■ As the project progresses there must be a concerted effort 

in establishing the criteria for the interaction, between the 
four physiological systems; respiratory , cardiovascular, ther- 
mal, and the renal-body fluid systems. Research such as that 

of Albergoni et al . (13) which attempts to define the signifi- 

\ 

r cant variables used m the interfacing model of the circulatory 
and respiratory systems and Lhe recently published research 
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effort of Guyton et al. (l4) will help in the formulation of 
these criteria. Guyton's research may prove most valuable 

i ___ * 

.since it attempts to integrate some of ,the dynamical inter- 
actions of various physiological subsystems; i.e, circulatory 

dynamics, pulmonary dynamics, regional blood flow control' and 

» 

its relationship to P , tissue fluids and pressures, heart 

°2 

rate, and stroke volume. 

To fully appreciate the transition between model choices 

1 

or, the potential capabilities of a model an understanding of 
the preceding references is mandatory. This is certainly 
not an exhaustive listing but several references have excel- 
lent bibliographies which makes it an adequate st at’e-of-art 

r 

listing for this research. 
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II. DESCRIPTION OF GR'ODINS ' RESPIRATOR! CONTROL MODEL 

The material presented in this section is a condensed 
summary of the model's description. Basically, this is the 
Version of the model that was used for the simulations involving 
the parameter sensitivity analyses with various exercise 
levels and environmental conditions. These simulations will 
■.be explained , ±n a later section. A more complete description 
of Grodins * model is included in the Appendix. 

Also, in.cluded in t.his section are comparisons of several- 
controlling functions and a brief statement of the immediate re- 
search goals as described ’in the following Study Report. 


A . General Model Description 

The model can be envisioned as a feedback control system 

comprised of a "plant 11 (the controlled system) and the regu- 
I 

latory component (controlling system). The controlled system 
is partitioned into three compartments corresponding to the 
lungs, brain, and tissue with a fluid interconnecting path 
representing the blood. A set of differential-difference 
equations describing a ga’s transport and exchange system 
between the lungs, blood, brain, and tissue compartments 
constitutes the framework for the models. A controller 


equation describing the inspired ventilation as a weighted 
function of the H concentration in the cerebrospinal fluid 



the H concentration m the arterial system 

t )), phe H concentration in the venous blood 
ao < 
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•of' the train (C ,„.+ i(t - x _ ) ) , and the oxygen partial pres- 

a \ n ) a a 

sure at the carotid chemorec eptor sites {F^q provides the 
feedback mechanism. A good qualitative explanation of the 
alveolar-arterial, venous blood-brain, and venous blood- 
tissue equilibria formulations involving acid-base buffering, 
metabolism, and empirical relationships is given by Grodins 
(8). These particular relationships which are used in the 
development of the differential-difference equations are pro- 
bably the most complex, both mathematically and physiologically 
speaking, of any of the model segments. The respiratory feed- 
hack system is represented *by the block diagram in Figure 1, 

, i 
i 

With reference to the' 1 , following description of the compart- 
• ment ali zat'i on and notation used, a list of symbols is given in 
,bhe Appendix of Grodins' paper (8). The same terminology is 
adhered to here and in the following discussions. The basic 
‘assumptions and features of the model's compartments that are 
characteristic of Grodins' program are enumerated below. 

Lung Compartment : 

1. Void of respiratory cyclic phenomena. 


2. Constant volume and uniform content. 


3. Flow rate of uni-form gas stream is dependent upon 
the variations ‘of RQ from unity, 

^a(C0 o }’ ^a(O-)’ an< ^ ^a(N_) of expired air and lung 
s c. d 

t 

exit arterial blood' are ' equal . 

Tissue Compar Linn n t. : • 

1. Uniform partial pressures equal to venous blood partial 


pressures . 



(8) 



C CSF ( H'*’ ) 

°a.(H + ) 

f a(o 2 .) 


Figure 1. 


Respiratory Feedback Control System 
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2. Removal of 'CO^ and supplying of 0 2 at a constant rate 
are functions of me.tabolism. 

3. Only physically dissolved 0^ and Ng are considered. 
^t(C 0 ) ^" S a ^' unc ^^ on °f (BECO^)^ and appropriate 
buffer relations. 

t . 

5. Bicarbonate’ concentration y H concentration. 

Brain -Compartment: 


1. Uniform partial pressures equal to venous blood 
partial pressures. 

2. Removing of C0 2 and supplying of 0^ at a constant 
rate are functions of metabolism. 

3. Only physically dissolved 0g and are considered. 

4. ^b(C 0 ) -*- s a function of (BHC0 3 ) b and appropriate 
buffer relations. 

5. Bicarbonate concentration ^ H + concentration. 

6. Brain/CSF membrane is permeable only to respiratory 
gases. 

* 

7. Partial pres.sure gradient occurs only across membrane. 

8. No buffering of carbonic acid in CSF , Bicarbonate 
content is constant for P CSF ( G0 ) > 10 mmHg. 

Blood Interconnecting Pathways : 

\ 

1. C /. is a function of P , n arterial. pH, and the 
a( 0g ; at 0 2 ) , 

Bohr effect (effect of ) ) . 


2 . G 


a(C0 ) is a function of P , (BHC0_) ,Hb, Hh0 o , 

d. Q* l LUg / 3 u c. 

plasnfa protein content, and the Haldane effect 


(effect of P 4 ( o 2 )). 



do) 


3. 


b. 


5. 


6 . 


T 


Arterial blood pH is a function of P . and bicar- 

a(C0 2 ) 

bonate content evolving from the Henderson-Hasselbaleh 
equation. 

Only physically dissolved nitrogen is considered. 
Transport delays are a function of volume and flow 
rates . 

Q and are functions of P , . and P , - . with 

B # a ( COg ) a(0g) 


assumed time constants. 

Recombination of venous tissue and venous brain blood 
enter the lungs after a variable transport delay. 


Controlling System ; 

1. H'ot explicitly described in terms of receptor and sensor, 
locations, afferent nerves, central nervous system, 
motor nerves, and respiratory muscles. 

2. Ho dynamically described elements in neurological 
aspects of the system — no time delays. 

- 3. Direct input (chemical concent rat ions ) --output (vent- 
ilation) relationships. 

1 < C 

An important aspect of the simulation is related to the 
blood transport time delays. As w-ith any physically dynamic 
system involving flow dynamics or transfer of information, there 
ar'e time delays in signal transmission. The time delay, t, 
is associated with the blood flow from the lung to the brain 

' 2 t 

compartment, from the lung to the tissue compartment, from the 
tissue to the lung compartment, from the. brain to the lung 
compartment, and from the lung compartment to the carotid 
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body receptor. These are very necessary terms and their incor- 
poration in the system equations enhances the confidence 
placed upon the dynamics of the simulation. To further refine 
these time delays, they are not considered to be predetermined 
constants but are continuously recalculated. They are time 
dependent in that they are functions of cardiac output, brain 
blood flow, and vascular segment volumes. This reasoning is 
,ver.y sound since the blood flow dynamics obviously affect 
relative site gas concentrations. However, consultation with 
Dr. Grodins has led to the questioning of the improvement in 
the simulation's fidelity with these time delay formulations 
.in comparison to> a more simplified version with a correspondingly 
simpler computer algorithm. 

A modification of Grodins' program involves the determin- 

t ♦ 

ation of respiratory frequency and heart rate. (15) At 
present the simulation is using the following expressions. 
Respiratory frequency is defined as 

FREQ = 8.1 d 7-815 (RMT(02) + RMB(02)) 
and hea'rt rate is 'defined as 

HRATE * 43.8 ( RMT ( 02 ) + RMB(02)) + 5^.5 
where HMT(02) and RMB(02) refer to the (^metabolic rates of 
the tissue and brain compartments respectively. 

B . Compan son o f Various Controlling Functions 

Tne controilmg function is a very important link in the 
total feeibach :r’.tam. This is one component of the study that 
has not been evaluated in depth other than variations of the 
parameter veightidgs j in the existing equation as i3 further 
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described in the parameter sensitivity analyses section. 

Other forms for the controller equation should he incorporated 
in Grodins ' model for a more complete evaluation. 

Since all controller functions simulate the actual system 
in a more realistic physiological manner over spe.cific variable 
ranges, it seems instructive to compare some of the more estab- 
lished ones . 


1 . Gray's Controller Equation 
Gray's controller , e quati on (7) describes alveolar ventil- 
ation as a function of alveolar CC» 2 , ^ACO ’ and a l ve °l ar 
oxygen, V AQ , as 


V A " aP A C 0 2 - b + C < d - P A0 2 > n i 0 

o 


( 1 ) 


with empirical constants, a, b, c, d, and n. is really a 

function of two expressions, one describing the dependency 
upon P AC0 , ( aP ACO “ b )> and the otber related to P AQ , 

, 2 jj 2 2 

(c(d - ) )• Consideration of only additive effects of 

C0 2 and 0^ has been proven invalid, but the inclusion of this 
equation illustrates one of the first stages of development of 

;a controller equation. 

i 

2 . Lloyd and Cunningham's Controller Equation 

1 

Lloyd and Cunningham's controller equation was developed 
under the basic assumptions that the product as well as the 
additive form of 00^ and 0„ effects be included, that exper- 
imental data be used. to empirically develop the equation 


for P 


ACO, 


'4 normal values, and that it would describe minute 
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ven'til itlon , V, It is given m tlie reference by Milhorn et 
al . ( 7 ) as 

? = D(P AC0 2 - B) [1 + A/(p ao 2 - C) 1 5 °- < 21 

The product of the two terms, D and the 'bracketed' expression, 

describe a hyperbolic expression with P. as the dependent 

variable * 

Here, P^qq 5. B = 37*24 rnmHg and P^ ^ C = 25 mmHg. 

1 ’In comparison to Gray's equation alveolar ventilation is given as- 

o' o o 

v A = v - (v/v T ) v D (3) 

with empirical expressions describing and as 

°b 

V .j, = V 1 = Tidal Volume 

and Vj, = Mg + bg = Dead Space, 

All values for the empirical constants are given by Milhorn 
1 (T)* A brief discussion describing the' salient features of 
Lloyd and Cunningham's equation and Gray's equation is also 
given by Milhorn (7). 

i 3. Milhorn' s Controller Equation 

i ■■■ 1 ** “*■ — ■' * ' ' '■ ~ rT ~‘ - • 

Although the controlling equation utilized in some of 

Milhorn' s more recent research appears to be of the form of 

Lloyd and Cunningham's equation, the reader is also encouraged 

to follow a more descriptive discussion of an earlier devel- 

opment in Milhorn's text (l6'V* The alveolar ventilation, V , 

is described as a function of the aortic-carotid bodies arterial 

Q_ concentration { C . . ) and the brain tissue C0 o concentration 

- - A ' iU o *=■ 

i— 

(Ctj ). It is given as 
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1 /. 


V A - a( °B 


CO, 


-b- + d(m-K 5 P B C A602 ) n , 0. 


(4) 


The term (C 


1 / 

) * 


C 02 


2 comes from an empirical CO^ dissociation 


curve. Also it is assumed that the relationship between 
arterial P_ and alveolar P„ is 

, i °2 2 

P a0 2 * "5 P A ' 


Then with 


K 5 F A 


— i c V C 
0 2 5 B Ag0 2 


(5) 


( 6 ) 


we have the form of the' last term in Equation 4. Other con- 
stants in Equation 4 are described as follows. First of all, 
i.t should be noted that some of the constants will assume differ- 
ent values for 1 the excitatory and inhibitory phases, 
a = 810.0 (Excitatory) 

a — 99.0 (inhibitory) 

b = 194.5 (Excitatory) 

b = 19.6 (Inhibitory) 

kg- 0.415 (Obtained from C0g dissociation curve) 
k^= 0.92 (normal arterial Pog/Normal alveolar Pq 2 ) 

-h-v 

. d = 8.0 (10) / Determined under inhibitory conditions 

n = 3.0 J 

\ 

m = 98.0 (Normal arterial Pq 2 ) 

4 . Grodins 1 Controller Equation 
The controller equation that is used in Grodins* program 
is given here with the same terminology as previously used and 
which ,is used in the Appendix. 
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v I = 1.138o[c(l^)C a(K + ) (t - ^ aB ) + (1.0 - C( 16 ) )c CSF(H +)] 

+ 1 *15^°C a ( H +) (t -< 0 ) + TERM - VI ( w ) (T) 

' where 

TERM = 23.6(l0)~ 9 [l04 - (B-47 )F a(q j(t “^o^" 9 

for (B - U7)P./ n At -< ) < 104 . (8) 

AlOgJ ao 

TERM = 0 for ( B 47 )F . , n , ( t - V ) * 104 . ( 9 ) 

, A(O^) ao 

The obvious question that needs to he answered is whether 
an expression like Equation 7 is more physiologically represent- 
ative than one like Equation 3. A summary of such a comparison 
yields the following list: 

1. Grodins’ modified program defines minute .ventilation as 
TVNT - DEADVT + (VE + Vl)/2. 

TVNT * ( .1107) (FREQ) + (.0?85)(VE) + ( VE + Vl)/2.(l0) 
Milhorn's use of Lloyd and Cunningham’s minute ventil- 
ation equation (Equation 2) should be compared to 
Equation 10. 

2. Grodins’ inspired ventilation, V^, corresponds to 
the alveolar ventilation expression used by Milhorn 

- which is Equation 3. 

3. One major difference between Grodins’ controller 
equation and Lloyd and Cunningham's controller equation 
is that wi&h the former there is no product effect 

of C0' o and 0 o . 

d. f. 

4. - Grodins' system his r,ho capability of distinguishing 

between the' sensing mechanisms of the CSF compartment 
and the 'vepous blood in the brain. 
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5. Grodins’ model reilec ts the metabolism and. acid-bise 
buffering by utilising a weighting of the H + concen- 
tration in the CSF compartment and at the carotid 
bodies site . 

6. Appropriate time delays are utilized for lung-to-brain 
and lung-to-carotid bodies in Grodins’ model. This 
feature simulates a more representative system in the 
dynamical sense. 

It is very difficult to justify the use of some 

of the -weighting factors in bo-th systems . For a more 

complete comparison, a form of Milhorn's controller 

equation should be incorporated in Grodins' program. 

-}* 

Then the benefits of H concentration sensing and/or 
blood gas tension sensing can be discussed m better 
perspective. This particular research will be 
elaborated upon in a future Study Report. 

C.. Re s e arc h Goals for the' Study Report 

As described in the contract's Statement of Work the 
underlying goal of this Study Report is to critically compare 
Grodins' respiratory program responses to Environmental 
Physiology Branch and published data and recommend and implement 
approved program changes to increase program fidelity. In 
order to fulfill v this goal a description of the mathematical 
development of Grodins' model is presented. See Appendix. 

.To evaluate and analyze the program responses a parameter 
sensitivity evaluation was performed using normal environmental 
conditions at a resting state and also environmental conditions 
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approximating those of Skylab with various levels of exercise. 

Since the controller equation was shown to be a very 

influential component of the system a special sensitivity 

+ , 

analysis for the weightings of the H concentration m the 
cerebrospinal fluid and blood compartments was formulated. 
These results are shown in the following report. When the 
exercise data from the Environmental Physiology Branch is 
received a more quantitative evaluation of the program's 
capabilities will be made. As it now appears the program 
responds with logical transient and steady-state trends for 
exercise levels. 
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I'll. VARIATIONS IN MODEL PARAMETERS 


A . General Considerations and Philosophy 

Before discussing specific parameter variations, a few- 
general ‘comments are appropriate. The primary overall goal 
of this investigation is to evaluate the model sensitivities 
for reasonable variations in parameter values. Grodins 1 mod- 
ified program .provides for various parameter variations by the 
use of input data cards. There are many standard conditions 
for inspired volumetric fractions of gases from which to base 
a sensitivity analysis. The decision was made to use the 
following data from which all simulations of variations in 
parameters are compared. 


B 

FI ( C02 ) 
FI (02) 
FI(N2) 


760 mmHg 
0.1000 
0.1100 . 
0.7900 


The time responses for V , Q, 


and 


P a(C0 2 ) 


are shown in Figures 


2-k respectively. 

It is believed that these base conditions are adequate for 
the sensitivity analyses performed. However, one should use 
extreme caution m projecting, m a direct fashion, the sensi- 
tivity of any parameter to other base conditions. Only one par- 
ameter or group of common parameters was altered at one time. 

f 

This is really the s first step in what should be a continuing 
effort in analyzing combinations^ of parameter 'variations and 
weightings . 


The format followed m the ensuing discussion is one which 
hopefully conveys the model's capabilities and limitations in a 
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concise and Informative manner. There is a tremendous amount 
of data which will not he discussed or presented in graphical 
form. This is not mean't to deemphasize their importance. Only 
the variables which illustrate significant changes will be addres- 
sed. The reader shouid acquaint himself with the output data 
format and all the possible relationships between variables that 
can be generated. See General Electric TIE 7^1-MED-3008 and 
TIE 7^1-MED-3009. 

With regard to inspired gaseous concentrations only the on- 
transient and related steady-state conditions are analyzed. 

The off-transient as related to the transition from exercise to 
the resting state is presented later in this section. Since 
the base run indicates an oscillatory response in the off-trans- 
ient of some variables, parameter sensitivity analyses will no 
doubt prove enlightning within this time range. 

The computer program is a modified version of Grodins ' 
.program similar to that used by Weissman (15). Although it has 
a subroutine which involves a work-exercise phenomenon, only the 
case of the resting subject is_ considered here. 

Table 1 presents the input data format giving card number, 
symbol, normal initial .value, and a brief discription of each 

- r 

parameter. In the discussion of each parameter variation,, the 
I 

card- number, symbol, and new initial value (or % change in base 
initial value) will be given. 



Card Mo. 

1 

2 

, 3 i 

4 
,5 
6 ' 

7 ' 

8 
9 

10 . 

11 


12 

13 

14 

15 

16 - 


17 : 


18 

19 


20 

21 
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TABLE 1: 

INPUT 

DATA CARDS 

Symbol Normal Initial 

Value 

Description 

F A ( C02 ) 

FA( 02 ) 
FA(N2) 

•i 

-.0527 

.1514 

.7959 


Alveolar gas fractions (dry), 
volumetric fraction of gas, 
dimensionless 

CB ( C02 ) 
CB ( 02 ) 
CB(N2) 

.6397 

.0011 

.0097 


Concentration of gas in brain, 
liters (STPD)/liter brain 

CT ( C02 } 
CT ( 02 ) 
CT(N2 ) 

.6132 

.0014 

.0097 


Concentration of gas in tissue 
compartment. Liters ( STPD)/liter 
tissue 

ft 

QB 

6.0000 

.7370 

[ 

Cardiac output blood flow, liters/ 
min. , cerebral blood flow, liters/ 
min . 

PCSF ('C02 ) 
PCSF ( 02 ) 
PCSF ( N2 ) 

47.8529 

36.0047 

567.4731 


Partial pressure of gas in 
cerebrospinal fluid compart- 
ment, mmHg. 

TMAX 

30.0000 


Length of computer run, min. 

CENT SENS PT 

0,0000 


Central Sensitivity Partition. 

Weighting of the H+ concentration 
in CSF with that of venous blood 
in the brain. With C(l6)=0, 
zero weight is given to venous 
blood at level of the brain and 
a weight of one is given to H + 
cone, in CSF. 

HB 

.2000 


Blood oxygen capacity, liters 
(STPD)/liter blood 

R1 

R2 

.1000 

.1000 


Time constants for cardiac output 
response (Rl) and cerebral blood 
flow response (R2) for changes 
'in blood chemical composition. 

CNT SENS COF 
CRTD BDY SCP 

1.1380 

1 . 154.0 


Controller sensitivity weight- 
ings, i . e . , 

1 

4 


V I - 1-1380 [0(l6)0 a(H + ) (t-x ii£ ) ♦ 


(i.o-c(i6))c csp(II +j3 


+ 1.1540 C ) +TERM- ’/ 1 ( N ) 

■ a { H ; ao 
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Card No. 

Symbol Normal Initial. Value 

Description 




1 

where 

t = Blood transport ’ delay from 

lung to carotid body, 

t = Blood transport delay from 
lung to brain , 

VI (N ) = C ( 37 ) 9 and 

TKRM = function of F,<. \. 

A( 0 g ) 

22 

KL 

3.0000 

Volumes of Lung (alveoli), brain. 

23 

KB 

' 1.0000 

and tissue compartments, liters 

21+ 

KT 

39.0000 


25 

MRB ( C 02 ) 

.0500 1 

Metabolic rates by brain. 

26 

MRB (02) 

.0500 

liters ( STPD ) / min 

2 T 

D ( C 02 ) 

;8i. 9900 

Diffusion coefficient for gas 

.28 

D ( 02 ) 

4.3610 

across "blood-brain", liters 

• 29 

D (N 2 ) 

2 . 524 o 

( STPD ) /min per mmHg 

30 

B 

760.0000 

Barometric pressure, mmHg 

31 

Fl( C 02 ) 

.1000 

Volumetric fraction of gas (dry 

32 

F I ( 02 ) 

.1100 

inspired), dimensionless. This 

33 < 

1 

FI ( N 2 ) 

. 7900 

is the value for the Base Run 
that’s utilized. 

34 

KCSF 

.1000 

Volume of cerebrospinal fluid, 
liters , 

35 

T 

.0000 

Init lal time . 

36 

f 

1 H 

.0078 . 

Size of computer time step, min. 

37 ' 

VI ( N ) 

87.5500 

Constant that is involved in the 
controller equation (See C( 2 l)).- 


Determines the normal level of 
alveolar ventilation so that 
?A( C 02 ) a ' S: rest, breathing 

air at sea level. When the con- 
troller sensitivity weightings 
are changed Vl(N) should be 
altered accordingly. 

38 Vl(SS) 5.3900 Value used for normal resting 

alveolar ventilation. This 
is not used in the program if 
VI ( N ) is known . 
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Card Wo . 

Symbol Normal Initial Value 

Description 

• 39 

PRINT AL TIM 

0.5000 

Output printed in these time 
increments. However, there 
is an over-riding statement 
that permits no increments 
greater than 0.5 min. 

40 

UNKNOWN 

1 

0.0000 

Importance related to C(39), 
but doesn't seem to be of any 
real significance. 

41 

BHC03 Blood 

.5470 

Standard bicarbonate' content. 

42 

BHC03 Brain 

.5850 

liters C0p (STPD)/liter x s 

43 

BHC03 Tissue 

.5850 

37°C where 

44 

BHC03 CSF 

.5850 , 

1 

t 

X - Blood, brain, tissue, CSF. 

45 

RMT ( C02 ) 

.1820 

Metabolic rates by tissue. 

46 

RMT (02 ) 

.2150 

liters (STPD)/min. 

4? 

DJI 

. ,0000 

Used in performing Dejours 

48 

DJ2 

.0000 

experiment (Not utilized in 
present runs). Brief description 



- 

of Dejours work relating 0^ and 




C0g threshold effects is given 
in Grodins 1 paper (8). 

'49 

*Col 1-6 W0RK2 

0 .00 

Work load (watts) 


Col 7-9 


Blank 


Col 10-15 DURAT** 




30,00 

Run time for work load 
(min. ) 


*F6.2 




**If DURAT is less than TMAX (Card 15) another work load 

* 

•» 

card is re'ad 

when print time 

exceeds DURAT. 

B. Con 

trolled System Parameter Var iatx ons-Barometric Pressure = ?60 mmHg 

Run 

No. Card No 

. Symbol 

Revised. Initial Values. 

107 C(l8) 

Rl 

. 0700 


C(19) 

R2 

.0700 

108 * c(l8) 

Rl 

.1300 ■ 


C(19) 

R2 

.1300 

Discussion: The values R1 and R? correspond to the time constants 

r l anci 

r n in the first 

order differential equations describing cardiac 

output 

and brain blood 

flow. By using 

70$ and 130$ of the normal initial 
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values essentially the time that is required for the blood flows 
to reach their steady-state values was altered. The model was 

shown to be insensitive to this alteration for t 5 . 1 min. There 

\ 

.were- slight variations in the initial transient (O^tsl), but 
even in this time period, the system seems fairly insensitive 
to these 30 % variations. Cardiac output reaches its steady- 
estate value in about 4 minutes, so in general its transient 
* period is less than for some ( other variables. 

Run No, Card No. ^ Symbol Revised Initial Values 


112 

C ( 4-2 ) 

BHC03 

BRAIN 

0 .5660 


• C(43) 

BHC 0 3 

TISSUE 

0.5660 

113 

■ C(42) 

BHC03 

BRAIN 

0.5470 


C ( 43 ) 

BHC03 

TISSUE 

0.5470 

114 

C(42) 

BEC03 

BRAIN 

0.5850 


C( 43 ) 

BHC03 

TISSUE 

0.5470 

115 

C(42) 

BHC03 

BRAIN 

0.5470 


c(43) 

BHC03 

TISSUE 

0.5850 


Discussion: Refer to Figures 5-7. Grodins indicates 

that the standard bicarbonate content of the brain and tissue 
compartments of O. 585 O is too high. Several simulations were 

I - „ 

run from which some data are plotted in Figures 5-7. Physio- 
'logically speaking there is a lower limit placed upon possible 
values for C(42) and c(43). Neither can be less than the value 
of C(4l), the bicarbonate content of the blood. This extreme 
limit was one of several chosen for evaluation. The model was 
'-■found to be much more sensitive to changes in the bicarbonate 
content of the tissues than of the brain. This is likely to 
be related to the relative volume sizes of the two compartments 
in addition to the buffering and metabolism relationships. The 
most pronounced changes occured in the transients with the 
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steady-state responses remaining unaltered except for the H 
concentrations related to the brain and tissue compartments, 
i«e« , and ^ ( j£ * 


Run No. 

Card No . 

Symbol 

Revised Initial Valu 

129 

C(27) 

D ( C02 ) 

50 # 

Normal = 40.9950 


C ( 28 ) 

C ( 02 ) 

50# 

Normal = 2.1805 

130 

C(27) 

D ( C02 ) 

150# 

Normal =122.9850 


C ( 28 ) 

D ( 02 ) 

150 # 

Normal = 6 . 5415 

131 

C(27) 

D ( C02 ) 

5 0# 

Normal = 40.9950 


C(28) 

D ( 02 ) 

50# 

Normal = 2.1805 


C(29) 

D(N2) 

50# 

Normal = 1.2620 

132 

C ( 2 7 ) 

D ( C02 ) 

150 # 

Normal =122.9850 


C ( 28 ) 

D ( 02 ) 

150 # 

Normal = 6.5415 


C(29) 

D(N2) 

150 # 

Normal = 3.7860 


Discussion: Refer to Figures 8-9. The qualitative trends 

that these variations produced are very physiologically sound. 

Noting that the diffusion rate is directly proportional to the 

diffusion gradient, we have a direct change in the rate of 

COg, 0g , and N^ passing across the brain-CSF barrier as we look 

at - 50 # and 150 # variations in their diffusion coefficients. 

Since the blood-cerebrospinal fluid barrier is very impermeable 

to H + , a change in the C has little effect upon the C„ c .. r ,/ TT +x 

^ a ^ n 1 bbf \ H ; 

However, p a ( C Qo) variations are detected in the CSF(l 8 ). The 

blood-cerebrospinal fluid barrier is highly permeable to 00^5 

thus, a 50# or 150# variation in D(C02) will be significant. 

In addition, when C0 o enters the CSF compartment, C_ oW „+x is 

£ U br \ n ) 

increased via the formation of carbonic acid. This is reflected 
in the Henderson-Hasselbalch equation and provides for the CSF 
compartment to be very sensitive to CO,, variations . One should 

u 

not expect a. significant change in for variations of D(02) 
or D'(N2) since M effects are not considered in the controller 












yr±li 
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equation and the simulated dependency upon 0 2 would not signif- 
icantly affect the ^C 3 p(j£ + ) or ^ a (H + ) directly controls 

Vj. The simulations support these conclusions in "being very 
sensitive to D(C02) variations and insensitive to D(02) and 
D(U2) variations. Both an increase in the transient and steady- 
state responses were obtained for and O CS p^ H +j with Run 
Nos. 130 and 132. Corresponding decreases were observed for 
Run Nos. 129 and 131. 

C. Controlling System Parameter Weightings-Barometric Pressure = 


760 mmHg. 


Run No. 

Card No . 

Symbol 


Revised Initial Values 

201 

C ( 2 0 ) 

CNT SENS 

C0F 

50$ Normal = 

0.5690 


c{ 37 ) 

VI ( N ) 



62 .6399 

202 

C ( 20 ) 

CNT SENS 

C0F 

90$ Normal - 

1.0242 


C( 37) 

VI (N) 



82.5680 

203 

c ( 20 ) 

CNT SENS 

C0F 

110$ Normal = 

1.2518 


c ( 37) 

VI ( N ) 



92 . 5321 

2 Oil 

C ( 20 ) 

CNT SENS 

COF 

150$ Normal = 

1.7070 


C ( 37) 

VI(N) 



112.4602 

Discussion: 

Refer to Figure 10 

. The variation in 

C ( 20 ) 


corresponds to a variation in the weighting of ^ CS p( H + ) in the 
controller equation. As previously described respiratory control 
is very sensitive to alterations in the C CS p^+^. There was a 
corresponding change made in C ( 3 T ) so as to compensate for main- 
taining P / n * & 40 mmHg at rest. If a new calculation for 

a \ t Ug / 

,Vl(N) was not made, the system would use the value of C ( 3.8 ) and 
then calculate a new value of Vl(N), This step was bypassed 
: with hopefully a good initial approximate calculation for Vl(N). 

■* No significant variations were observed in any compartmental 

■f 

H concentration, Q, or p a ^QQ )• However, as illustrated 















(35) 


in Figure 10, V^ is directly dependent upon C(20) although not 
responding rapidly in the initial transient. 


Run No. 

Card No . 

Symbol 


Revised Initial Values 

205 

C ( 21 ) 

CRTD BDY 

SCF 

50$ 

Normal 

* 0.5770 

206 ‘ 

C(37) 

vi(n) 




65.9842 

C ( 21 ) 

CRTD BDY 

• SCF 

90 $ 

Normal 

= 1.0386 


C ( 37 ) 

VI (n ) 




83 .2368 

207 

C ( 21 ) 

CRTD BDY 

SCF 

110$ 

Normal 

= 1.2694 


C( 37 ) 

VI'(N) 




91.8632 

208 

C ( 2 1 ) 

CRTD BDY 

SCF 

150$ 

Normal 

= 1.7310 


■ C ( 37 ) 

VI ( N ) 




109.1158 


Discussion: Refer to Figures 11-12. With consideration 

given to all the system variables *the simulation seems to be 
more sensitive to variations in the carotid body sensitivity 
■weighting, C(2l), than the central sensitivity weighting, C(20). 
Here, the initial phase of the transient is modified with this 
modification significantly realized in the cardiac output, 

0<t$4 min. Although the steady-state values of the other var- 
iables are not altered appreciably, the system's initial trans- 
ient variations in V and Q are reflected in the other variables' 
transient responses. The term Vl(N) is used in the same manner 
as in Run Nos. 201-204, Physiologically speaking, one should 
expect the system to respond faster to changes at the carotid 
receptor sites than at the CSF site. Not illustrated in such a 
pronounced fashion, but certainly evident in Run Nos. 205-208, 
is the fact that <the system that has the faster .response seems 
to reach a steady-state value of lower magnitude than does the 
system with the slower response for all variables except V^ 


and Q. 
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n No . 

Card No. 

Symbol 


Revised Initial Values 

209 

C ( 20 ) 

CNT SENS 

COF 

50$ 

Normal = 0 . 5690 


C (21) 

CRTD BDY 

SCF 

50$ 

Normal = O .5770 


C(37) 

VI ( N ) 



41.0741 

210 

C ( 20 ) 

CNT SENS 

COF 

50$ 

Normal = 0 . 5690 


C( 21 ) • ‘ 

CRTD BDY 

SCF 

150 $ 

Normal - 1.7310 


C(37) 

VI (N) 



84.2057 

211 

C ( 20 ) 

CNT SENS 

COF 

150$ 

Normal - I. 7 O 7 O 


C( 21 ) 

CRTD BDY 

SCF 

50$ 

Normal = 0 .5770 


C( 37) 

VI (N) 



9Q .8944 

212 

C ( 20 ) 

CNT SENS 

COF 

150 $ 

Normal = 1.7070 


C ( 21 ) 

CRTD BDY 

SCF 

150 $ 

Normal - 1.7310 


0(37) 

VI ( N ) 



134.0260 


Discussion: Refer to Figures 13-15. There are some inter- 

esting interacting phenomena that occur -when Doth C(20) and 
C(2l) are varied. From initial investigation of Figure 13, 
one might conclude that for both C(20) and C(2l) equal to 50)5 
or 150 $ of their normal values that the system’s responses are 
additive. However, when the system's responses for independent 
variations in C(20) and C(2l) are viewed and compared to Figure 
13, it is noted that the responses are not additive, hut convey 
interaction. A decrease in C(20) by 50$ is more influencial 
than a concurrent increase in C(21) by 50$. In a similar manner, 

'an increase in C(20) by 50$ is. more dominating than a concurrent 
decrease in C(2l) by 50$. As steady-state conditions are approached 
the simulation becomes more dependent upon the CSF detecting 
system than the carotid body site. The opposite effect can 

I 

be mentioned for the initial part of the transient. Here the 
carotid body site is more influential. 

Figure 14 shows the effects of the extreme conditions upon 
cardiac output. Figure J 5 showing the response of ) bo 

the extreme variations of Ull'O) nud 0(2 L) i lino bruboa that the 
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faster responding system eventually reaches a lower steady-state 
value than the slower responding system. The variations in 


3 a(CO ) eady-stat e 


steady-state P 


a(C0 2 ) 


are in agreement in that a larger 
maintained when V ^ is reduced and 


vie e- 


ver sa . 


Run Ho . Card No . Symbol Revised Initial Values 

116 C ( 16 ) CENT SENS PT 0.1000 

117 C(l 6 ) . CENT SENS PT 0.2000 

118 C(l 6 ) CENT SENS PT 0.4000 

119 C( 16 ) CENT SENS PT 0.6000 

120 C( l 6 ) CENT SENS PT 0.8000 

Discussion: Refer to Figure 16 . The central sensitivity 

partition term is involved in a weighting of the response to 
H concentrations in two compartments, namely the CSF and train. 
Although CSF H concentration is the dominant controller in the 

train-CSF compartment, it has teen suggested that venous train 

+ 

H concentrations are also involved in the controller function 
for ventilation. In addition, the latter responds more rapidly 
than the CSF detector system. With C(l 6 ) = 0, zero weighting 
is given to the venous tlood at the level of the train and a 

m Hb 

weighting of one is given to the H concentration in the CSF. 

The form of the controller equation reflecting this condition 
is shown as 

■ vi = i.i380(c(l6)c a(H+) (t.T aB ) ♦ ( 1 . 0 - 0 ( 16 )) c CSF(H+) ) 

+ 1. 15*10 0 ) + TEEM - VI(H). 

q» v xi ) a,o 

All of the variables illustrated a faster response as 
C(l 6 ) was increased; however, a variation in steady-state re- 
sponse was realized only in the inspired ventilation. This 
feature is a direct reflection on the fact that in the steady- 
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■f m t 

state response the H concentration in the arterial compartment, 

C /, T +\ = 59*5 nanomoles per liter Hood, is less than In the 

a(H )ss 

CSF compartment, ^cSF(H + )ss ~ ?1 • 5 nanomoles per liter CSF, 

D . Comments on Simulated Skylab Conditions 

The importance and usefulness of the model will be more 
evident in environmental simulations approaching the conditions 
of Skylah . From the simulations run for sea level conditions, 
barometric pressure = 7^0 mmHg, it was concluded that the most 
sensitive component of the system was the controlling equation. 
Thus the efforts of this section are concentrated on the 
controlling equation and its variable weightings both for the 
resting state and various exercise levels. 

The first simulation effort involved a change of only the 
barometric pressure and inspired gaseous concentrations. Thus, 
these specific input data cards were altered to read as follows 
with the COg concentration corresponding to 5 Torr, 

Card No . Symbol Revised Initial Values 


C(30) 
C ( 31 ) 
C ( 32 ) 

C ( 33) 


B 

FI ( C02 ) 
F I ( 02 ) 
FI(N2) 


260.0000 

.0192 

.7000 

.2808 


With the above conditions for the simulations, the subject 
was introduced into the new environmental conditions with no 
preconditioning period. Consequently, no faithful results 
could be associated with the transient solution. The steady- 
state conditions that were obtained for all the variables were 
of value since they were required to establish the appropriate 
initial values for all the simulations using the above conditions . 
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•"In a sense, the use of these steady-state values compensated 
for the preconditioning period. The transient period for most 
of the variables appeared to be less than 5 minutes; therefore, 
a simulation run of 25-30 minutes produced good steady-state 
values. Table 2 shows the comparison of the initial values 
for the variables at normal environmental sea level conditions 
versus Skylab conditions. The variables associated with'N^ were 
'altered the most. 

TABLE 2: COMPARISON OF INPUT DATA CARDS FOR SEA 

LEVEL VERSUS SKYLAB CONDITIONS 

Card No. Symbol Initial Value for Normal Initial Value for 

Environmental Conditions Skylab Conditions 


1 

FA ( C02 ) 

.0527 

.1783 

2 

FA ( 02 ) 

.1514 

.5336 

3 

FA(N2) 

.7959 

.2881 

4 

CB( C02 ) 

.6397 

.6413 

5 

CB ( 02 ) 

. 0011 

.0012 

6 

CB(N2) 

.0097 

.0011 

7 

CT ( C02 ) 

.6132 

.6153 

8 

CT ( 02 ) 

.001b 

.0015 

9 

CT(N2) 

.0097 

.0012 

10 

Q 

6.0000 

6.0000 

11 

QB 

.7370 

.7496 

12 

PCSF( C02) 

4? . 8529 

48.1202 

13 

PCSF ( 02 ) 

36.0047 

36 .6316 

Ik 

PCSF ( N2 ) 

567.4731 

70.6804 

15 

TMAX 

30.0000 

30.0000 
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Initial Value for Normal Initial Value for 
Card No. Symbol Environmental Conditions Skylab Conditions 


16 

CENT SENS PT 

0 . 0000 

. 0000 

IT 

HB 

.2000 

.2000 

18 

R1 

.1000 

.1000 

19 

R2 

. J 000 

.1000 

20 

CNT SENS COF 

1.1380 

1.1380 

21 

CRTD BDY SCF 

1.1540 

1.1540 

22 

KL 

3.0000 

3.0000 

23 

KB 

1.0000 

1.0000 

24 

KT 

39.0000 

39.0000 

25 

MRB ( C02 ) 

.0500 

.0500 

26 

MR-B ( 0 2 ) 

.0500 

.0500 

27 

D ( C02 ) 

81.9900 

81.9900 

28 

D ( .0 2 ) 

4.3610 

4.3610 

29 

D(N2) 

2.5240 

2.5240 

30 

B 

760 . 0000 

260.0000 

31 

FI ( C02 ) 

.1000 

.0192 

32 

FI (02) 

.1100 

.7000 

33 

FI(N2) 

.7900 

.2808 

34 

KCSF 

.1000 

.1000 

35 

T 

.0000 

.0000 

36 

H 

.0078 

.0078 

37 

VI ( N ) 

87.5500 

87.5500 

38 

Vl(ss) 

5.3900 

5.3900 

39 

PRINT AL TIM 

.5000 

. 5000 

40 

UNKNOWN 

.0000 

.0000 

4l 

BHC03 BLOOD 

. 5470 

.5470 

42 

BHC03 BRAIN 

. 5850 

.5850 
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Initial Value for Normal Initial Value for 
Card No . Symbol EnTxronniental Conditions Skylab Conditions 


43 

BHC03 TISSUE 

• 5850 

. 5850 

4 4 

BHC03 CSF 

.5850 

. 5850 

45 

RMT ( C 02 ) 

.1820 

.1820 

46 

RMT ( 02 ) 

.2150 

.2150 

4? 

DJI 

.0000 

.0000 

48 

DJ2 

.0000 

.0000 

One 

of the goals of the 

Skylab experiments 

is to evaluate 


the subject's performance at various exercise levels. The 
base run to which the succeeding parameter sensitivity analyses 
were compared incorporates the exercise phenomenon. Exercise 
levels of five-minute intervals were used. The work loads 
and corresponding time intervals were described as 

Time Interval Work Load * 

0 <t< 5 min 0 watts 

5 min $t< 10 min 50 watts 

10 min <t< 15 min 100 watts 

15 min £t< 20 min 150 watts 

20 min $t< 30 mm . 0 watts 

At 5} 10, and 15 minutes an on-transient of all the variables 
may be observed. The of f-transient or return to the resting 
state is contained in the interval from 20 to 30 minutes, 

E . Controlling System Parameter Weightings - Barometric 
Pres sure = 260 mmiig 

Although significant results were obtained for all system 
variables only a representative number of outputs were plotted. 
Inspired ventilation, cardiac output, an example of c ompar tmental 
* Note card No. 49 on page 25 for input data format 
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H concentration ( C„ „„ , „+ v ) , and an example of C0„ partial 

liSiU ) 2 

pressure in the blood (P , >) were used. Also data were 

a ( C Og J 

plotted illustrating the interrelationship between Og partial 


pressure m the blood (P 


Run No . Card Ho . 


303 

304 


C(16) 

c(l6) 


, „ v) and cardiac output, Q . 

a l Ug ) 

Symbol 


CENT SENS PT 
CENT SENS PT 


Revised Initial Value 
0.4000 
0.8000 


Discussion: Refer to Figures 17-20. From Figure 17 one 

can see that when more weighting is given to the H + concentration 
in the venous blood of the brain than in the CSF compartment 
there is an increase in oscillation in the Vj. response to 
exercises. All exercise levels produce a slightly lower steady- 
state response for when CENT SENS PT is increased. The 
decrease in exercise level to the resting state is not signif- 
icantly altered by variations in CENT SEN'S PT . Cardiac 
output. Figure l8, is unaffected except possibly for the off- 
transient where a slight oscillation is observed. 

The oscillations are more pronounced for 20-$ts30 minutes 
in Figures 19-20. Throughout the entire simulation an increase 
in the CENT SENS PT correspondingly raises the H concentration 
in the CSF compartment and the partial pressure of CO^ in the 
"blood regardless of the exercise level. 

In going from the resting state to the first exercise 
level the arterial C0 o tension is reduced indicating that the 
increase in inspired ventilation for this level of exercise is 
more than adequate. However, as the exercise level is increased 
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beyond, the 50 watts both the arterial C0g> tension and H + concen- 
tration in the CSF compartment is increased. This Justifies 
the statement that the ventilation rate is not as efficient 


for increased exercise. Also from Figure 17 it is evident 
that the ventilation rate increase- does not respond as swiftly 
for higher exercise levels. The more weighting that is given 
to the H + concentration sensing in. the venous blood of the brain, 
the higher the relative arterial CO^ tension is for any instant 
of time. As expected the level of the H + concentration in 
the CSF compartment is higher when the CENT SENS PT is increased. 


Run No . 

Card No. 

Symbol 

Revised Initial Value 

313 

C ( 20 ) 

CNT SENS COF 

50% Normal = O. 569 O 


C(37) 

VI ( N ) 

62.6399 

314 

C (20) 

CNT SENS COF 

150# Normal = 1.7070 


C ( 37) 

VI ( N ) 

112,4602 


Discussion: Refer to Figures 21-24. There was no signifi- 

cant deviation from the base run for the variations in the 
weighting of the H + concentration in the CSF compartment. A 
slight variation occurred in the arterial COg tension response 
when the initial exercise level was implemented. Under the 
existing environmental conditions (low inspired C0^ concentration) 
the system is insensitive to various weightings of ^csf(H + )* 

Of course this fact is evident since from Figure 24 it is 


shown that c C s F (h + ) not 


deviate very much from the base 


run . 


pi f J/v5.^ 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Run No . 

Card No. 

Symbol 

Revised Initial Value 

323 

C ( 2 1 ) 

CRTD BDY SCF 

! 

50$ Normal = 

'“ 0.5770 


C(37) 

VI ( N ) 


65.9842 

324 

C ( 21 ) 

CRTD BDY SCF 

150$ Normal = 

1.7310 


C ( 37 ) ' 

VI ( N ) 


109.1158 


Discussion: Refer to Figures 25-28. The arterial C0 2 

tension m Figure 27 displayed the most sensitivity toward var- 
iations in the weighting of the H + concentration at the carotid 
bodies. An increase in the CRTD BDY SCF weighting produced 
a more pronounced oscillation in the transient response of V^. 
when the fixst exercise level was simulated. Basically, the 


c-arotid bodies' receptor site is a faster responding system 
component than the CSF component; therefore, it is not too 
surprising to observe slightly more sensitivity for similar 
variations in CRTD BDY SCF than in CNT SENS COF. 

Many times it is more instructive to analyze the relation- 
ship between two variables instead of viewing the variables 

on a time basis. The two variables chosen are P \ and Q 

a \ Og / 

as shown in Figures 29 and 30. The two curves can be explained 
by following through a sequence of events corresponding to the 
letters A through H. Figure 29 corresponds to a deemphasis 
in the H + concentration weighting at the carotid bodies' 
site while Figure 30 corresponds to an increase in the same 
weighting factor. 

Time is measured along the curve connecting the letters 
with t=0 at point A ana c=3Q minutes at point H. The sequence 
of events which relate to the exercise changes can be described 


as follows. 
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Interval A-C : 5~'1Q' minutes. Work load = 50 watts 

, 4 * 

Interval C-E: 10-15 minutes. Work load = 10Q watts 

Interval E-G: 15-20 minutes. Work load » 150 watts 

Interval G-H : 20-30 minutes. Work load - 0 watts 

(return to resting state). 

Intervals A-B, C-D, and E-F correspond to very short 

time segments. All three of these intervals are less than 

\ minute in length. As the exercise level increases, the 

cardiac output increases and the partial pressure of. 'Og 

In the hlood decreases. Both Figures 29 and 30 illustrate 

a recovery in P , n % for, the exercise levels. The steady- 
a(0 2 ) ■ 

state values of P , „ % and Q for the three exercise levels 

a (, 0 g / 

correspond to points C,E, and G. When the maximum exercise 
level is reduced to the resting state, the two variables take 
on values along the curve from point G to point H. 

The recovery phenomenon referred to in the previous 
paragraph can be physiologically described as follows. Basically, 
an increase in exercise promotes an increase in cardiac out- 
put. The increase in cardiac output along with an increase 
in the ventilation rate is not adequate to maintain the arter- 


ial tension at a constant level for the exercise steps of 
100 and 200 watts. The increase in ventilation rate is more than 


ad eq.ua t e for the first exercise leve 1 s i-nce there is an increase 
in both cardiac output and arterial 0g tension. After the car- 
diac output approaches its steady-state value for each exercise 
level the arterial 0 o tension recovers as shown by the curves 
approaching points E and G. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 

Upon performing some preliminary studies as described in 
the, previous sections, it appears that Grodins * Respiratory 
Control Model can be made adaptable to the research efforts 
of the Environmental Physiology Branch of the Biomedical Research 
Division, NA.SA-MSC. There are several facets of the program 
that need to be altered which will provide for a more complete 
physiological representation. On the surface there appears to 
be potential benefits from the incorporation of other research 
efforts, i.e., Milhorn’s models, while utilizing the basic 
framework of Grodins* model. This phase of the project will be 
considered later. 

I 

The parameter sensitivity analyses need to be extended to 
include combinations of parameter variations and also to criti- 
cally evaluate the of f-trans lent (t 5 30 min) conditions. A 
test for the physiological constraints of the model with regard 
t'O C0 2 input levels would be an important investigative effort. 
Levels of inspired CO^ are physiologically correlated with 
yentilation rates, coma, anesthetized states, and death. The 
model should' perform in a similar manner. If not, COg limits 
must be specified. 

The most flexible feature of the entire system is the con- 
troller equation. There needs to be a continuing effort in the 
modification of the existing controller equation and in the eval- 
uation of other controller equations. This is especially im- 
portant in view of integration of the respiratory system model 
with other physiological systems. The controller equation will 
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be a very important link in the total system providing some of 
the interface capabilities between the cardiovascular, thermal, 
and body fluid systems with the respiratory system for both the 
resting and. exercising subject. 

To perform with any acceptable degree of fidelity, the 
system must be closely correlated to the salient physiological 
functionings. Variations in the solubility coefficients of the 

* i , *r 1 

gases for each compartment and a reduction of the bicarbonate 
.content of the brain and tissue compartments , seem to be desir- 
able. The model' includes no general description of the control 
of cardiac output and regional blood flow. Certainly the regional 
,hlood flow- feature will be' important when the system is inter- 
faced with the body fluid and cardiovascular systems. 

As mentioned in Section II. A, respiratory frequency which 
is now defined as 

FREQ =.8.1 + 7.815 (RMT(02) + RMB(02)) 
should be expanded to illustrate other physiological dependencies 
in addition to 0 2 metabolic rates of the tissue and brain com- 
partments. In a similar manner, heart rate, defined as 
HRATE = U3.8 ( RMT ( 02 ) + RMB(02)) + 5^.5 
should be expanded to include relationships involving stroke 
volume, regional blood flow, increased work load, and perhaps 
a proper weighting of parasympathetic and/or sympathetic stim- 
ulations . 

Probably the most important element of the simulation is 
the adaptation of the model to particular individuals. This 
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will require a better defined set of initial parameter values 
in addition to relating these parameter values to the desired 
environmental conditions. 

At present there is a feeling of skepticism with, regard 
to the exercise portion of the simulation in that it doesn't 
relate to the actual physiological processes involved. There 
is supportive experimental evidence that indicates the initial 
transient ventilation response due to a change in the exercise 
level under stable environmental conditions is created primar- 
ily by neurological control (l8). With the present simulation, 
a change in the exercise level causes a change in the metabolic 
rates of CO^ and 0^ in. the brain compartment which via chem- 
ical control is reflected in the ventilation rate. A control 
mechanism which appropriately weights the neural and chemical 
processes is being studied. This postulated change will also 
be related to the functions describing heart and respiratory 
rates during exercise. 

Investigation will continue in the forementioned areas 
with the goal being to develop a model that is physiologically 
sound and compatible with the research of the Environmental 
Physiology Branch. Specific system alterations will be given 
upon further evaluation of experimental data related to 
exercise levels and after a more thorough investigation of 


Milhorn’s research. 
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V. APPENDIX 

The following material is a description of the mathematical 
development of Grodins ' program (8) excluding the exercise 
simulation component. The description involves 

(A) Compartment material balance equations 

(B) Alveolar-arterial concentration equilibria 

(C) Venous blood-brain equilibria 

(D) Venous blood-tissue equilibria 

(E) Cardiac output and brain blood flow 

(F ) Blood transport time delays 

(G) Controller equation 

(H) Differential-difference equations. 

All of the explanations including equation numbers refer to 
Grodins’ paper (8). 

A . Compartment Material Balance Equations 

Equations 1.1-1. 3 of Grodins' model (8) will not be 
duplicated. This format will be carried throughout the ensuing 
discussion. These three equations describe the time rate of 
change in the nondimens ional volumetric fraction of COg , 0^ , 
and in the lung compartment. They are first order differ- 
ential equations containing the following terms. Each equation 
contains the difference between two product terms. One of the 
terms is the product of the inspiratory gas flow rate multiplied 
by the volumetric fraction of each gas in dry inspired air. 

The other term is similar in that it is the product of the 
expiratory gas flow rate multiplied by the volumetric fraction 
of each gas in dry alveolar gas. Evident from these two terms 
is the fact that water vapor is not considered in this gaseous 
exchange re lat i on snip . The difference between the venous and 
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arterial concentrations of each gas is multiplied by the\ 
total cardiac output.- One should note that the concentration 
terms have units of. liters per liter bLood. Another coefficient 
relating barometric pressure and water vapor pressure of the 
blood multiplies the concentration 'differences, All of the 
terms are divided by the volume of the .lung <?ompart.ffl&ht 
yielding the proper volumetric fraction of -gas per minute. 

I 

Expressions of the time rate of change of each gas in 
the brain and tissue compartments are defined by Equations 
1.4 - 1.9. The metabolic rates of CO^ and 0 2 have opposite 
effects in both the brain and tissue compartments. Both are 
considered constants in each compartment. The metabolic 
rate of CO increases the rate of change in concentration of 
CO^ while the metabolic rate of 0^ decreases the rate of change 
in concentration of 0^. There is no metabolic rate term 
associated with the nitrogen gas. This concept corresponds 
to the basic assumptions in the previous sections. The differ- 
ence between the arterial and venous concentrations of each 
gas in the brain compartment, Equations 1.4 - 1.6, is multi- 
plied by the blood flow in the brain, thus altering the rate 
of change in the concentrations of each gas. Also, within 
the brain compartment there is a membrane separating the, CSF 
subcompartment from the brain. The membrane is permeable 
to all three gases, being about 20 times more premeable to 
CO., than 0 o as depicted by the' relative magnitudes of the 

k. «— 

diffusion coefficients. The product of the difference 
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between the partial pressures of each gas in the brain and 
CSF and the diffusion coefficient correspondingly alters 
the rate of change in concentration of that gas in the brain 
compartment. The brain — CSF interface is described in the 
manner that contributes to a positive alteration in the rate 
of change in the concentration of a gas if the partial pressure 
of that gas is greater in the CSF compartment than in the 
brain compartment. 

In addition to the metabolic term described above, the 
rate of change; in concentration of each gas in the tissue 
compartment. Equations 1.7 - 1-9, is a function of the tissue 
blood flow. ' Tissue blood flow is defined as the difference 
between cardiac output and brain blood flow. This difference 
is then multiplied by the gas concentration difference between 
the arterial and venous sides of the tissue compartment. Once 
again, the volumes of each compartment are entered into Equa- 
tions 1.4 - 1.9 to provide the proper dimensions of liter 
per liters of blood per minute. 

The partial pressure of each gas in the CSF compartment 
can be described by a first order differential equation as 
in Equations 1.10 - 1.12. The time rate of change in the 
partial pressure' of each gas is a function of the diffusion, 
coefficient across the blood-train barrier, the volume of - 
the CSF compartment, a conversion factor relating atmospheres 
and ramllg, the solubility coefficient of the gas in the CSF, 
and the difference between the partial pressure of the gas 


across the blood-brum barrier. 


As expected, the dimensions 
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of Equations 1.10 - 1.12 are mmKg per minute.. It should "be 
noted that the solubility coefficient of each gas is given 
in units of liter per liters per atmos'phere and evaluated 

O t - 

at 37 C. These coefficients are not altered between compart- 
ments and satisfy Henry f s Law for the environmental conditions 
being considered. The constant k is a conversion factor. 

From the -initial assumptions associated with the lung 
.compartment it is obvious that a change in the volumetric 

1( t 

fraction of one gas is related to the sum of the changes in . 

, the' volumetric fractions of the other two gases. Therefore, 


^'he- algebraic sum of + F A ( q 2 ) + F 

leads to Equa v txton- 2.1(8).. Both V £ and 


a(n 2 ) 

are 


is aero , which 
needed for the 


solution of the system. Thus, Equation 2.1 gives a relation- 


ship from, which V ' can be written in terms of V T or vice versa, 
v I* 1 ** , ht x 

, , - » 

Vj is an output of the controlling system. The expired 
ventilation, V , is a function of the concentration differences 
of each gas multiplied by the total cardiac output and a 
non-dimensional coefficient related to the barometric pressure. 

This is a reasonable relationship if one considers that 
the output flow rate is equal to the input flow rate plus 
alterations that might occur in the metabolic processes re- 
lating” C0 o , 0_ , and . The latter is tempered by the cardiac 
,>, 22.2 

ou-tput and the difference between the input-output coneentra— ; 
tions of each gas., 
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B . Alveolar - Arterial Concentration Equilibria 

There are probably several methods for building up the 
necessary interrelationships between variables so as to solve 
the final system of differential-difference equations. Grodins ' 
next step in the development is tp formulate expressions for 
the alveolar-arterial concentration equilibria for the gases 
in the blood at the lung exit (Equation 3.1 - 3.7). The 
Haldane effect which illustrates the dependence of the con- 
centration of COg upon variations in Og concentration is 
included. Also included is the description of the dual role 
that variations in C0 o concentration play upon the concen- 
tration of Cg. This latter dependency is termed the Bohr 
effect . 

The cone entration of COg, Equation 3.1, is comprised of 
a constant bicarbonate term, the last term being a weighting 
of the effect of the partial pressure of COg, and the two 
middle terms being chemically combined relationships. Although 
not readily explainable because of the coefficients involved, 
the second term is a weighting of the difference between the 
blood oxygen capacity and the concentration of oxyhemoglobin 
(liters of 0 per liter blood). The third term attempts to 
describe the product effect of the presence of Og and COg with 
the Og associated with the Hb term and COg described by the 
logarithm of CQg arterial concentration and dry alveolar 
gaa concentration difference divided by a weighted arterial 
partial pressure of CO;,. 
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The expression for C & ( Q Equation 3.2, contains an empirical 
relationship describing the 0 o dissociation curve and. the Bohr 
effect. The concentration of oxyhemoglobin is described by 
an empirical form containing an exponential term which is 
dependent upon the pressure of CC> 2 whose involvement becomes 
evident through the pH term. An explanation of the derivation 
for pH is in order. The normal form for the Henderson- 

a. 

Ha^selbalch equation is 


pH = pK + log / HC0 3 


Grodins has defined pH as 

a 

pH = - log C 
* a s a ( H ) 

with the units of being nanomoles per liter. Thus 

pH a = - log [C a(H+) do)' 9 ] 

pH a = - log [C a(H+) ] . log ri0 - 9] 

P H a 9 - log [ c a ( H + )]- 

To relate this development to the Henderson-Hasselbalch equa- 
tion requires the terms of Equation 3.6 to be regarded as physi- 
ological equivalent to the terms in the Henderson-Hasselbalch 
equation , i ,e . , 


pK <=> log K 

and HC °3 < ~ > ^ a C0g (B-4 t) Eft(COo) 

C ° 2 °a(CO a ) - Ho tB-I.T)P A(c 0 , 

c- C 


These analogies have been developed by several researchers for 
various body compartments; one example being given by Bradley 
et al . (19 ) . 
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Since only - 'physically dissolved nitrogen is considered, 
the arterial concentration of N 2 . can be defined by the 

simple expression in Equation 3.7. The development describing 
the, alveolar-arterial concentration .equilibria ana the develop- 
ment in the following section concerning the blood-brain equil- 
ibria are probably the most complex, both mathematically and 

physiologically speaking, of any segment of the respiratory 
control model. 

G * Venous Blood-Brain Equilibria 

The expressions of the concentration of C0 2 in the brain 

and Tenous blood reflect the original assumptions in that the 
brain concentration is not dependent upon Hb and HbO £ effects,. 
While the venous blood is described m exactly the same manner 
.as the alveolar-arterial CC> 2 concentration. One other variation 
. 1S that the stan dard bicarbonate content differs in the two 
^expressions (Equations U.l.atid 4.2). The bicarbonate content 
of the brain, (BHCO^, .must always be greater than or equal 
to the bicarbonate content of the blood, (BHCO^. This partic- 
ular relationship between these two parameters is one which ■ 

• heeds to be .altered so as to more precisely simulate physiolog"- 
ical conditions. 

The 0 2 concentration in the brain, C b(q is given as 
C B(0 2 3 = ka B(0 2 ) P B(0 o )‘ 

This agrees with the form for C B( } vhteh is given as 

^ Cb ( w 2 ) = U b(h 2 ) P b(w 2 )- 

These t*t> relationships are utilised because in the brain 
ment only physically dissolved 0 g and W 2 are considered. 


compar t 
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The concentration of 0 ^ in the blood at the brain exit is 
also a function of the oxyhemoglobin concentration. This parti- 
cular involvement is described by Equations ^.3-4.7. As in 
the previous section, the oxyhemoglobin concentration is an 
empirical form containing an exponential term which is dependent 
upon the venous blood pH, pH . As before, to relate this de- 
velopment to the Henderson -Hasselbaloh equation, from Equation 
fc>7 


and HCO. 


<=> 


CO, 


pK <*=> log K ' 

kct co 2 p b(co„) 


vB(C0 2 ) - ka co ^ P B ( C0 ^) 


The first expression in Equation i.8 states that the concentration 
of in the venous blood leaving the brain is proportional to 
the concentration of in the brain. However, in the present 
simulations the coefficient is unity since the solubility coefficient 
any gas is the same for all compartments. This is another rela- 
tionship between parameters that needs to be altered to more 
clbsely represent the physiology involved. 

D . Venous Blood - Tissue Equilibria 

As is suggested by Grodins (8) the CO^ , 0 ^ , and relation- 
ships between the tissue compartment and the tissue venous blood 
are derived under the same basic assumptions as those of the 
venous blood-bram equilibria. Thus one can develop an analogous 
set of relationships for the tissue compartment as for the 
brain compartment. Since the physically dissolved concepts, 
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the metabolic, relationships, and the empiricisms are analogous, 
then the overall system sensitivities to variations within 
the two compartments ■ must be closely tied to i the relative 
magnitudes of t'he c ompartmental volumes. This conjecture appears 
to he justified' in some of the computer , simulations . 

E . Cardiac .Output and Brain Blood Flow 

Both cardiac output and brain blood flow are defined to 
be dependent upon variations in the partial pressures of arterial 
‘Og and COg . The first order differential equations that des- 
cribe these two flow patterns are similar. Looking at the 
differential equation for cardiac output. Equation it 

t 

can be rewritten as 


- r i « + « - V 

Th!e function, f n , can "be considered the 1 forcing function for " 

r * ^ 

the differential equation. Note that by Equation 7*2, f is a 
, * Q 

i 

collection of terms dependent upon the normal cardiac output, 

a variation in the blood flow due to Og and a corresponding 
5 1 . 

variation due to COg . Let's ignore any physiological signifi- 
cance of .the dit f erential equation. The solution, Q(t), can 

be considered as the sum of two components, Q (t) and Q (t).' 

c P 

Q ft ) is the solution of the differential equation if f =Q. 
c « Q 

.Such a solution would take the form of 
’’ Q c (t) = Ke " t /r 1 

where K vould’ be determined from the initial .value Q(t), 

The other component of the solution, Q (t), is determined from 

. P 

the' form of the' forcing function f^. A simple example would be 
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to ignore any effects of 0^ and CO^ and consider f = . 

, > » 

By the classical solution of differential -equations, it is 
readily shown that 

Q (t) = Constant * Q„. 

P N 

The total solution would he 

Q(t ) ■= Q ( t ) + Q (t) 
c P 

or Q(t) - Ke -t /ri + Q w 

Suppose that at time t = 0 seconds, Q(t) = Q(o). Then 

t=0 

K = Q ( 0 ) - Q m . 

and the total solution would he 

Q(t) = ( Q ( 0 ) - .'“ t /ri + 

A graphical description showing two values of Q(o) relative to 



Figure 1. Sample of first order response 
Up to this point the mathematical description of cardiac 

-output depends upon the value of the initial condition, the 

* 

normal resting value, and a time constant which attempts to 
drive the solution to a physiologically meaningful representation. 

f 

.Fortunately, hy the addition of an empirical relationship 



both defined over 


and a linear relationship describing ^ 

specified limits for P |l(o } and * a(co ) 


re-spec tively , Cfrodins has mathematically formulated a more 
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precise description of Q(t}', There is a factor AQ ^ ^ contri- 
buting to the forcing function only when' j <10^ mmHg. 

1 , > 2 

The computer program deviates from the paper in the math- 
ematical description of the dependency of AQ^^ ^ upon P a (QQ ) 
Equation 7.6 “has been altered to read 

AQ( C0 j = 6 for Y 5 60 

and AQ (C0 2 ) = 0 for Y S 40 . 

These factors are described by Equations 7.3 and 7.5. 

The same form of differential equation describes the brain 

blood flow except that the forcing function components related 

to P , nr . \ and P \ are different. The same range of P /_ ■* 
aiCOgj aiOgj akOgj 

is considered, but the empirical relationship is different 

than the one used for cardiac output. On the other hand, 

P , . n is a contributing factor for all of its possible vsj-ues 
a l o ’ • 

except for 38 mmHg g p a (co ) * ^ In this range, brain 

blood flow is not dependent upon P ». Variations in the 

a i u Og ) 

time constants r^ and r^ are considered in some computer simu- 
lations involving parameter sensitivities. 

F . Blood Transport Time Delays 

As with any physically dynamic system involving flow 
dynamics or transfer of information, there are time delays m 
signal transmission. The time delays of the respiratory model 
are' associated with the blood flow between compartments. The 
notation used for the time delays can be readily described by 
Figure 2 which indicates the shape of the delayed waveform and 
the notation used. 



(80) 



In the simulation, the time .delay, t, corresponds to delays 
between lung and brain compartments , between lung and tissue 
compartments, between tissue and lung compartments, and between 
the lung compartment and the carotid body receptor. These are 

very necessary terms and their incorporation in the system 

> 

equation enhances the confidence placed upon the dynamics of 
thq simulation. Concentrations of the three gases at the 
entrances of the brain, tissue, and lung compartments are defined 
in terms, of the appropriate time delays by Equations 8,1 - 8.9. 

To further refine these time delays, they are not considered to 
be predetermined constants , but are time dependent in that they 
ard functions of cardiac output, brain blood flow, and vascular 
segment volumes. This reasoning is very sound since the blood 
flow dynamics obviously affect relative site gas concentrations. 
These time delays which are continuously being recalculated 
in the computer program are defined by Equations 8.10 - 8.l4. 

From the integral representations it is clear that when steady- 
state conditions are reached, i.e., Q and Qg become constants, 
the time delays will .become, constants also. This is also 
justified by observing the computer print-out. 

G. Controller Equation 

[i — — 1 — — ‘ ' 

The controller equation is the system's feature that pro- 
vides th'e feedback mechanism. Many 'different forms for the 
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controller equation can be postulated. To completely under- 
stand the simulation's potentials, a concentrated effort 
should be given to investigation of other controller equations. 
This reasoning is amplified by Grodins ' ^comments and conse- 
quently special emphasis is placed upon its evaluation in the 
study. At this point one might ash the following question. 

What is the purpose of the controller equation? This question 
might be answered best by observing the block diagram of a respir- 
atory feedback system in Figure 3. 



F A(0 o ) 

C. 


Figure 3. Respiratory Feedback Control System 
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point the thesis is that, C cg p ^ , c a ( H + p and ^a(0 ) 

» , 2 


The controller equation as described by Equation 9*2 uses 
only three of many outputs of the controlled system. At this 

are the ■ 

three factors which determine the Inspired ventilation, V^. . 

The inspired ventilation is then used to determine and 

is also used in the controlled system's equations. Thus, the 

feedback path is completed. Also, it is very instructive to 

note that the arterial H + concentration and alveolar oxygen 

volumetric fraction are not detected until t minutes after 

ao 

■the arterial blood leaves the lungs. This feature simulates 
the time elapsed while the blood is flowing to the carotid 

4 > 

receptors from the lung compartment. The H concentration in 
the CSF subcompartment is described by Equation 6.1. In the 
same manner as the pH calculations are performed in the other 


compartments a comparison of terms in the Header s on-Has selbal'ch 
equation with those of Equation 6.1 yields 

pK <=> log K* 

and HC0 3 ^ ka CSF(.C0 2 ) P CSF(C0 2 ) 


CO, 


(BHC0 3 } CSF 


Once the decision has .been made as to the variables which 

> i 

are used as inputs to the controlling equation, a proper weighting 
of these variables is made so as to drive the total system's 
responses to correlate with experimental data. This necessarily 
requires enroirical formulati on$ . Although not explained in the 
paper the computer program, nas the capability of weighting the 
influence of the H+ concentration in the CSF compartment to the 



(83) 


H-+ concentration of the venous blood of the brain compartment. 
This weighting is evaluated in some of the computer simulations. 
Great care should be expended so as not to critically violate 
known physiological features. 

H « Di f ferenti al-D i f f erence Equations 

It is now apparent that many mathematical relationships 
have been developed relating alveolar-arterial concentration 
equilibria, venous blood-brain and venous blood-tissue equil- 
ibria, blood flow rates to all the compartments, and the trans- 
port time delays. The next task is to assemble these relation- 
ships in a meaningful set of equations from which solutions 
for the system's variables can be obtained. Grodins (8) 
hag accomplished this by utilizing the framework of the compartj- 
mental mass balance equations. Although it is not necessary to 
completely describe the culminated effort in Equations 10.1 - 
10.9 and Equation 11.1, perhaps it would be instructive to 
comment upon the substitutions for one of these equations, i.e., 
Equation 10.5. Equation 10.5 has the basic structure of Equa- 
tion 1-5. For convenience these two equations are repeated here. 
Equation 1.5: 


B <° 2 ) " K 


B 


-MR 


B 


(0 p ) + Q B( C aB(0 ? ) " C vB(0 )) ~ D 0 p ( P B ( 0 ) “ P 


OSP 


Equation 10.5: 

o 

b( 0 2^ ic 


B 


" MR B(0 2 ) +Q 3 


r fa 02 \ 

~(aiTo 2 )J 


fK 

_ui 


[A] 

i 


2 


C B(0 o ) ~ C vB(Hb0 


2>) ( 


p b(o 2 )~ p csf(o 2 )J 



(8b) 



is substituted,. This substitution corresponds to bracketed 
term [b]. All other terms remain the same in forming— the 
differential-difference equation . 
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INTRODUCTION 


I . 

A group of respiratory system models which is of interest 
from the viewpoint of evaluation for their application within 
the goals of the NASA project is that of Milhorn and models 
related to Milhorn's research. In a previous Study Report 
a modified version of Grodins 1 model was investigated illus- 
trating its response to particular parameter variations and 
various work load conditions. Host of Milhorn' s models, 
although developed under different experimental conditions 
and for varying motivations, contain features that are similar 
to Grodins 1 work:. For this reason it seems appropriate to 
discuss the major features of Milhorn' s models and perhaps 
suggest special modifications for Grodins' model which might 
utilize these features. Therefore, the major thrust of the 
following material will be to summarize the theory and 
assumptions of Milhorn 1 s models keeping in mind the conditions 
for which they are applicable and to describe potential 
applications of Milhorn 's work to improve the existing Grodins' 


program . 



II. MODEL VARIATIONS AND SIGNIFICANT FEATURES 


Within the listed references ( 1-6 ) Milhorn et al. have 
published several model variations each of which was developed 
for a particular experimental arrangement or environmental 
condition. Instead of correlating the results of the referenced 
research and confusing the various system terminologies it 
will be more instructive to independently deal with each 
research effort. This approach will be more justifiable as 
the constraints and experimental developments are enumerated. 

In addition, the models will be discussed in a chronological 
manner. 

A. Milhorn' s Basic Compartment al Respiratory Model 
Many of the assumptions and mathematical formulations 
of the compartment al model of the respiratory system that 
appears in Milhorn' s text (l) are also used in his more recent 
work. Here it was assumed that the chemical factors that 
control alveolar ventilation are CO^ and 0^ with the effects 
of the H+ stimulus related to P by a linear relationship. 

The basic assumptions can be summarized as follows (l): 

1. The circulatory system containing no mixing phenomena 
but with regionally associated time constants inter- 
connects the lung, brain, and body compartments. 

2. Dissociation curves for 0 o and CO. are equal for 

C. d. 

arterial and venous blood. 

3. Brain blood flow is dependent upon cerebral-arterial 

P„_ and relationships. 

co 2 o 2 
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it. Arterial P n is linearly related to alveolar P_ . 

U 2 U 2 

5. Venous F is equal to tissue P~ . 

°2 U 2 

6. Arterial P nn is equal to alveolar P nn • 

7* Venous P is equal to tissue P nn • 

C0 2 t0 2 

8. Linear relationship between P(-, A and P which 

L0 2 °2 

controls minute alveolar ventilation. 

9. Respiratory quotient equals unity. 

10. Metabolism and cardiac output are constant. 

11. System simulates acid-base balance disturbances. 

12. H+ concentration is a function of P co . 

13. Carotid and aortic bodies are sensitive to 0^ . 

The model comprizes a controlled system and a controlling 
system. External disturbances in either the form of CO^ 
excess, 0^ deficiency, or both are considered inputs into the 
controlled system. The regulation aspects of the_sy_stem 
are accomplished via ventilatory and circulatory parameters. 


P P 

a_ * a n , and H+ are monitored outputs of the controlled 
C0 2 °2 

segments of the system. These variables are weighted giving 

o o 

a value for alveolar ventilation ( V ) . V^ is then considered 
a second input to the controlled system along with the 
environmental disturbances, thus completing the negative 
feedback formulation. 

Input-output block diagrams are utilized in the development 
of the equations which quantitatively describe the system. 

There is one such diagram for each of the following components . 
1. Ventilatory controJler 


I 


2 . 


Lung compartment 



k 


3. Brain compartment 

U. Body compartment 

5. Venous mixing 

6. Cerebral circulation 

7. Comparator for cardiac output, cerebral blood flow, 
and body tissue blood flow 

8. Transportation lags in circulatory system 

9. Overall compartmental model. 

The expressions for the time rate of change of COg 
and Og volumes in the brain, body, and lung compartments are 
developed. These six equations contain twenty-eight variables. 
Using the previously listed assumptions, mathematical relation- 
ships are derived such that eight variables are used in eight 
differential equations describing the normal human respiratory 
control system. The eight variables are 


1. 

q b 

O 

2. 

V A 

3. 

c B 

B C0, 

If. 

C T 

^CO 

5. 

C A 

co 

6. 

\ 

7. 

V 

8. 

°2 

C A 

°2 

Six 

of t 


cerebral blood flow (liters/min) 
minute alveolar ventilation (liters/min) 
brain tissue COg concentration (vol fraction) 
body tissue CO,-, concentration (vol fraction) 
alveolar COg concentration (vol fraction) 
brain tissue Og concentration (vol fraction) 
body tissue Og concentration (vol fraction) 
alveolar Og concentration (vol fraction). 


chan 


ge in the variables , C , C , C, , C , and C. 

r\ X _ X) X 


CO, 


'CO, 


CO, 
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Another equation is called the cerebral circulatory controller 
equation and describes the time rate of change in brain blood 
flow. The final equation, the ventilatory controller equation, 
describes the alveolar ventilation in terms of various 
weightings on the controlling functions of the medullary 
respiratory center, carotid and aortic bodies, and the 
mechanical portion of the lungs. In effect, this latter 
equation, described as 1_ 

v, = a(C ) 2 -b +d (m - rf. P-0. a ) >, 0 (l) 

A b co 2 5 E A fo 2 

provides the simulation with the necessary feedback mechanism. 

As was observed with Grodins * simulation the ventilation 

controller equation is a very sensitive feature of the entire 

system. Ventilation is controlled by a different set of 

variables in Milhorn’s model than in Grodins 1 model. However, 

the major difference evolves from the fact that Milhorn’s 

model simulates the ventilation as a sum of the effects of 

P nA , H+ , and P with the final programmed equation having 
LQ 2 °2 

the effects of the H+ concentration embedded in the empirical 
constants of the equation. To evaluate and compare the 
physiological responses of Grodins* simulation for a ven- 
tilatory equation like Milhorn’s equation and Grodins' 
equation one must be restricted to the resting state and 
normal sea level conditions. Increase in ventilation due 
to exercise would be determined only by the alteration in 
metabolic x’ufces and would not involve any neurological control 
as does the exercise controller equation of the present Grodins* 


model . 
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Sufficient experimental data was not available at the 
time; therefore, a description of the cerebral blood flow 
in Milhorn's model was based upon several assumptions. 

The first assumption is that a three-fold increase in oxygen 
concentration in air doesn't alter the cerebral blood flow. 

A decrease in oxygen concentration to 10 % increases the 
cerebral blood flow by 35$. For a further decrease to 5 % 

C >2 was assumed to double the normal cerebral blood flow. An 
empirical relationship describing the effects of P n 
upon cerebral blood flow (developed by Kety and Schmidt (j)) 
was utilized. Combining the effects of CO^ , O^, and normal 
cerebral blood flow, the equation is written with the 
empirical constants as 



Considerable emphasis was placed upon the steady-state 

analysis of Milhorn's model. The only transient analyses 

performed with model were on- and off- transient responses 

for alveolar ventilation as the system was subjected to 

instantaneous changes in the levels of C0 2 and 0^ concentrations 

at the inputs to the lung compartment. Probably the most 

beneficial result of the f or einent i one d research was the 

realization that much more experimental and simulation research 

was needed to fu Lly describe tin? respiratory control system 

and realise a simulation ' potentials. This includes the 

investigation of cerebrospinal fluid influence, variations in 

cardiac output and circulation time delays , and the relationships 

between alveolar ana end- alveolar capillary P_ and P PA • 

t u_ 
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B. A Model Illustrating Static, Dynamic, 

Mechanical, and Cyclic Phenomena 

In a progress report concerning a conceptual view of 

a digital computer model of the human respiratory system (2) 

Milhorn discusses "both the static and dynamic features of the 

model and builds upon the previous work. Within the static 

component the controlled system equations describe the alveolar 

tension of C0 Q and 0^ . Both tensions are described as functions 

of their respective inspired gas tensions , production and 

consumption rates of C0 2 and 0^ respectively, the respiratory 

quotient, and the barometric and water vapor pressure. The 

two controller equations considered were those of Gray, 

O m 

V. - aP. -1 + e (d - P. ) >, 0, (3) 

°°2 \ 

and Lloyd and Cunningham, 

V = D (P ~B) I 1 A / (P -C) I ; P } 40 mmHg. (U) 

A C0 2 L a o 2 J a co 2 

As described by Milhorn, 

V A = ; - fV D (5) 

where fV^ corresponds to minute dead space. All parameters 
described in this paragraph can be altered and appropriate 
responses observed. 

The dynamics of this model centered ground a restructuring 
of the ventilatory controller equation. Initially, experi- 
mental evidence indicated that ventilation response could 
be simulated by weighting the influence of the H+ at 55/» and 
the influence of C0 2 at This weighting scheme was 

further refined by experiments from Lambertsen 1 s laboratory 
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which associated special weightings of the H+ contribution 
from the peripheral sensors, a central receptor in the 
brain, and a central receptor in the cerebrospinal fluid. 


The equation illustrating this weighting is 

V = 0.12 V(C ) + 0.44 V(n ) + 0.44 V(C nQ _ 

PS H+ B H+ CSF H+ 


). ( 6 ) 


The contribution from the 0 2 sensing mechanism evolved from 
Equation 4. 

Exercise alters the gaseous exchange mechanism in the 
lungs. The effective area of the pulmonary membrane is 
increased by the change in blood flow. Milhorn has investi- 


gated this area change hy utilizing the equation. 


A = A °+ 1 f k-. VL + k_ Q 

k i + M \ r 

' ° C 0 


having the empirical constants, k and k 0 , V being the 

-L id Li 

o 

lung volume, Q c the pulmonary blood flow, and the subscript 
0 referring to normal values. No evaluation statements 


were made with regard to the simulation of Equation J. 

Milhorn suggests an analog flow diagram representing 
the mechanical dynamics of the system. No simulation results 
were available for evaluation of this model. Although the 
simulation was not detailed for the exercise phenomena a 
cyclic model was proposed in order to describe the possibility 
of the peripheral chemoreceptor s being unidirectional rate 
sensitive. It was lelt that the tie between the unidirectional 
rate sensitive detectors and the fluctuating arterial 

v 

blood H+ , CC > 2 , and 0 2 concentrations could be better defined. 
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However, other sources do not emphasize these fluctuations; 
therefore, this refinement might not he appropriate with 
relationship to other system components. In any event, the 
cyclic pulmonary system is described by'varying passageway 
resistances, anatomical and effective deadspaces, a-v shunts 
for the blood, and the feature of venous admixture. 

In this reference Milhorn also emphasizes the importance 
of the controller equation. The interrelationships between 
the cardiovascular and respiratory systems were experimentally 
investigated with the goal being to expand the range of 
Lloyd and Cunningham's controller equation range for P. 

co 2 

The experimental apparatus of Holloman, et al. (8) was 

utilized to establish an empirical model for P. < 40 mm Hg . 

C0 2 

Further results for P < Ho mm Hg was to be presented 

A C°2 

at a later date . 

C. A Model Emphasizing Venous Admixture 
Milhorn and Pulley ( 3 ) have published a theoretical 
study which is believed to be of value in the evaluation of 
diseased states of the respiratory system, i.e. emphysema, 
pulmonary congestion, or any condition that would affect 
’the pulmonary membrane's functions. One feature of this 
model is its description of the capillary shunts and the 
corresponding description of the blood gas tensions (venous 
admixture case) at the mixing poinL of the end alveolar 

the art-ovi'il blood (through the bypass). 


capillary blood and 
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OF POOR QUALm 


This study was not intended to he a complete simulation 
of the respiratory system, hut instead, a component was 
developed which would be incorporated m a more complete 
computer simulation. The development of the pulmonary gas 
exchange equations evolved from a description of the rate of 
flow of gas across the surface area of n parallel cylinders 
each representing a capillary segment. Here, the total 
number of capillaries is defined as 

n = (Q - Q )/f (8) 

S 

o 

where Q = cardiac output (lit/min) 

O 

Q g = pulmonary shunt flow (lit/min) 
and f = blood flow through a single capillary (lit/min). 
Further development yielded an expression for capillary gas 
exchange as a function of gas tensions. Capillary gradients 
were then developed dependent upon surface area, membrane 
thickness, cardiac output, or all possible combinations of 
these variables. 

As stated previously, one of the salient features of 
this simulation was that of the venous admixture phenomenon. 
Under normal venous and alveolar gas tensions it was shown 
that approximately 1.3 per cent of the total cardiac output 
was shunted past the capillary system. From this fact the 
simulation of the admixture phenomenon can be detailed. 

Although the simulation is adaptable to the study of 
inert gases, i.e. N ,, and lie* , the immediate concern is 
related co extrac taol>' !«aia v <-s aosociuted with 0 and C0 o , 
Interrelationships beiween C0 o and 0 o , the Haldane and Bohr 

i— d 
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effects, complicate the development. Considering the compart- 
mentalization of Grodins' model and the mathematical development 
associated with the capillary gas exchange mechanism of the 
model discussed here it is doubtful that any of this model's 
features can he directly applied to Grodins' system. The 
features of this model would he of value if a refinement of 
Grodins' model included a more detailed description of the 
gas exchange mechanism associated with the lung compartment. 
Another aspect of the model by Milhorn and Pulley (3) is 
that it illustrates no simple adaptation to the study of 
exercise and its control of ventilation. The most important 
advantage in the use of this model would he the component 
of the simulation associated with the admixture phenomenon 
as the fractional gaseous concentrations of the environment 
were altered. 

D. A Model Maximizing Steady-State Information 
In the steady-state simulation of the human respiratory 
system by Milhorn and Brown (k) emphasis is directed toward 
justifying that ventilation he described by a function of the 
receptor sites' stimuli as CO^ and 0^ levels are altered. 

This is extremely important if farthful reproduction of 
transient responses are to be obtained. Since receptor 
site locations, numbers ol sites, and the general nature of 
the responses at these sites are ill-defined this research 
emphasizes the steady-state respiratory system and attempts 
to more fully uii-lize t. lie Lugo quantities of information 
from the steady-state. 
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The system is divided into a controlled system component 

and a controlling system component. As described in some 

of Milhorn’s previous -work (l) the controlled system has 

o 

alveolar ventilation, , as an input calculated from the' 
controlled system vith alveolar 0^ tension and alveolar C0 2 
tension as outputs. These two outputs are treated as inputs 
to the controlling system, thus completing the feedback 
system . 

This particular model by Milhorn is more compatible 
with Grodins ’ model than the one previously described (l) 
in that there are system descriptions that might be readily 
adaptable to alterations in Grodins' model. Alveolar 0 and CQ^ 
tensions are described by two equations which contain variables 
which are closely associated with those of Grodins' program. 
These equations are repeated here, described in terms of 


inspired 0^, 0 consumption, alveolar ventilation, respiratory 


quotient, water vapor pressure, and barometric pressure. 


> = p _ V 

ao 2 io 2 o 2 


P H 0 + P I0 

d d 


3 , - P + V 

aco 2 I00 2 ’co 2 


P B ' P H,0 + PjCO, * 11 (10) 

d. d I\ 


Plots ol these equations are given for various values 


of Pjq and PjQ n With 2.5 £ < 100 liters /min and constant 

2 2 

0 n consumption. When P_ p is varied P. is low for all 

° o 2 

values of P^q if is low. As is increased a 

steady-state value of P, in assumed which approximates P Tn 

* 0 ., I0 2 
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for V A > ^ liters/min. The opposite relationships are 

o 

established between V , P , and P 

A it '° 2 A C°2 

When Equations 9 and 10 are plotted for various levels 

o o 

of 0g consumption (V ) with 2.5 5 < 100 liters/min, 

P Sttd F 

°2 ' C0 2 

for all levels of 0g consumption, if V A is low P^ is low and 

O 0g 

P. is high. As is increased beyond approximately 


vary in the following manner. As expected, 

. O • 


CO. 


6o liters/min a steady-state value for P (l^Q mm Hg) 

^2 

is assumed irregardless of the level of 0^ consumption. 


Likewise, as is increased the value of P^ 


decreases 


CO, 


rapidly and assumes steady-state values below 10 mm Hg 


for V A > 6 0 liter's /min. 


upon the 


The graphical dependence of P. and P. 

o 2 a co 2 

barometric pressure, P^, is shown for steady- state values 


of V., 2.5 S V S 100 liters/ min. P ana P 

a o 2 a co 2 

by a family of curves for variations in P„ 

±5 


are described 


(300 j Pg $ 1500 mm Hg) which are similar to the family of 

O 

curves ‘for V n variations. For these plots P Tn and P Tnr . 

°2 * I0 2 2 
were held constant at 149-3 and 0.3 mm Hg , respectively. 

O 

The relati onshins for P and P. with V. for different 

a o 2 a co 2 

values of the respiratory quotient (R) are plotted. For 


low values of V., P is low and P is high independent of 

°2 C0 2 

the specific value of*" R. However, as V is increased, 

A 

P is increased and approaches a steady-state value of 
°2 

1^0 mm Hg for >. h 0 liters/min. In a similar , 

O 

manner, for all values- of R , P is decreased as V increases 

A C0„ A 
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Although perhaps difficult to justify quantitatively, the 

O 

■above results for P and P versus V. as the different 

°2 ' Ac0 2 

variables in Equations 9 and 10 are independently 

altered are physiologically justifiable in a qualitative sense. 

The major thrust of Milhorn and Brown’s work with the 
controlling component of the model is in the comparison of 
Gray’s classical equation and Lloyd and Cunningham's equation 
(4)'. Gray's controller equation is very similar to Equation 1. 

t 

With the empirical constants evaluated it is written as 


o 

V, 


r if. 9 

= 1. 8p - 6 7.^5 + 4.24(10)"° ( 104 - P, ) z0 (ll) 


CO, 


Equation 11 is developed using assumptions that are not 
just i f lab le . The fact that CO^ and 0^ effects are not just 

additive as illustrated m Equation 11 has been shown. 

Also, experimental data was not available for the entire 
range of ventilation even though the equation was developed 

O 

to simulate the entire range of V^. 

Describing CO^ and 0 2 effects upon ventilation as 
multiplicative as well as additive is one of the major 
contributions of Lloyd and Cunningham's work. Experimental 
data shows that the relationship between ventilation and 
P can, be described by a straight line with a slope dependent 

co 2 

upon P . This slope was then described as a hyperbolic 
°2 

function of P, . Thus the interdependence of these 

°2 

variables are more involved than in an additive arrangement. 

The equation that evolved with values for the empirical constants 


is given by 
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v = 2(P -37.2)4) (1 + 13 . 6 / ( P . - 25)) * o 

• c °a . \ 


( 12 ). 


for P * 31.2k, 

• C0 2 


; The relationship between V and is developed as follows. 


where 

and 


V A = V " f V D 
= dead space 


(13)' 


v/v n 


respiratory frequency 


Dead .space is calculated from an empirical relationship. 
v D = o.o?85V t + 0.1107 (14) 

where V^, = 0.134 V = tidal volume. 

In addition to the equation that describes the relation- 

1 

ship between alveolar ventilation and minute volume other 
important qualities which distinguish Lloyd and Cunningham's 
equation from Gray's equation are listed here. 

• 1. As mentioned previously, CO^ and 0^ effects are 
both multiplicative and additive. 

2. It is not extended beyond known experimental data 

(valid only for P. z 37*24 mm Hg). 

C0 2 

3. Multiplicative effects produce a more significant 
response under the conditions of low 0^ - high CO^ 
levels . 

Several three-dimensional plots are presented by Milhorn 
and Brown illustrating the comparison of Gray's controller 
equation and Lloyd and Cunningham's control equation. The 
following combinations of variables were considered. 
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Case 

A. 

p io 2 » 

P IC° 2 - 

o 

' V A 

Case 

B. 

p i° 2 - 

p 

IC0 2 ’ 

V 

Case 

C. 

p 

I°2 » 

p 

ico 2 ’ 

f 

Case 

D. 

P A ’ 

P A 

O 

V A 



°a 

C0 2 


Case 

E. 

P A • 

s 

P A > 

a co 2 

V T 

Case 

F. 

p 

r A » 
°2 

P A » 

co 2 

f 


A very brief explanation of each case is instructive. 

o 

Case A: Peak values for occurred for low and high 

Pj^Q with Lloyd and Cunningham’s equation giving 

the larger values. Lloyd and Cunningham's results 

were undefined when P Tnn and- P-r^ were both low. 

ICOg IOg 

Case B: Peak values for occurred for low Pjq and -high 

■^ICO Lloyd and Cunn-ingham 1 s equation giving 

the larger values, Lloyd and Cunningham's results 

were undefined when P T _ and P Tn were both low. 

. Xu Ug -*■ ^2 

Case C: The highest frequency ( bre aths /min ) occurred for 


Case D 


low P and high P Tr , . Lloyd and Cunningham's 
results were undefined when P^ and Pj C q were both 
low . 

0 \ 

Peak values for V occurred for low P. and high 

°2 

P. with Lloyd and Cunningham's equation giving 

a co 2 

the larger values. Lloyd and Cunningham's 


results illustrate an undefined region for 
37.2*+ mm Hg . 


CO, 
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Case E: 


Very similar results 

equations for P. i 

C0 2 

Cunningham ' s 


for V,^ in the extreme 


? . and high P 
A i 


CO, 


were obtained for V^, 
37.2^ nun Hg. Lloyd 
equation did give 1 
cases corresponding 


with both 
and 

arger values 
to low 


Case F: Very similar results were obtained for f with both 

equations for P. 5 37.21 mm Hg. Once again Lloyd 

co 2 

and Cunningham's results were undefined for 

P . < 37-2lmmHg. 

C0 2 

A suggestion for implementation of Lloyd and Cunningham's 
equation into Grodins * model will be presented in a later 
section. Of course, this would have to be expanded to include 
the exercise phenomenon since this equation doesn't involve 
all of the physiological variables associated with exercise. 


E. A Model Evaluating Transient Ventilatory Responses 

Reynolds, et al. (5) developed an experimental technique 

for evaluating the transient ventilatory response to various 

mixtures of CO^ in air. The recorded outputs for the system 

included respiratory frequency (f), tidal volume (V ), 

alveolar oxygen tension (P ), and alveolar carbon dioxide 

°2 

tension (P ). Minute ventilation xn liters/min 

o C0 2 

(V) was " computed as a function of tidal volume and 

frequency . 

Because of insufficient technology to adequately measure 
the cyclic phenemona associated with ventilatory transients 
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at the time of these experiments, the time constants associated 

with the respiratory variables were not well documented. Also 

the correlation between the strength of the stimuli and the 

magnitudes ■ of the responses were uncharted. Therefore, the 

basic purpose of these experiments was to accurately record 

f , V , P , P , and V for different step inputs of C0 p - 
°2 C0 2 

air mixtures in order to better define the transient 

response and justify the respiratory control system models in 
existence . 

The experimental arrangement will not be completely 
discussed here. It should be mentioned that the experiments 
consisted of a 20-minute prestimulus period, a 25-minute 
stimulus period, and followed by an of f-tr ansient period of 
20 minutes. The stimuli were 3, 5, 6, or 7% C0 2 - air mixtures 

A trigger circuit detected the flow signal. The gases 
were continuously sampled and analyzed by a Beckman LB-1 
infrared C0 2 analyzer and a Thermox 0 2 analyzer. The 
expiratory flow signals, 0^ tension, CO^ tension, and the 
output of the trigger circuit were all recorded on magnetic 
tape. After the experiments the results were analyzed 
yielding the respiratory frequency, minute ventilation, and_ 
alveolar concentrations. 

A summary of the results is directed tovu.'i the categories 
of tidal volume, respiratory frequency, minute /entilation, 
alveolar CO , alvcol ir 0 transient half-times, and the 

d ci 

subjective effects for the subjects. 
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1. Tidal volume: Rapid initial transient -with a more rapid 

-off-transient with a'll stimuli. No overshoot in either 

on- or off-transient. 

2. Respiratory frequency: Slower on- and off-transient 

than with tidal volume. These data can he used in 
defining a better frequency evaluation in Grodins * model. 

3. Minute ventilation: Smooth response for both on- and off- 

transient with no overshoots. Steady-state values are 
reached in a shorter period of time than that for the 
steady-state values of frequency. The opposite is true 

in comparison to the steady-state values of tidal volume, 
h . Alveolar COg:' Very rapid on-transi ent s and off-transients 
with an overshoot occurring in the off-transient for 
all stimuli cases. This is in agreement with simulations 
using Grodins* program. Lower CO^ levels also caused 
overshoots in the on-transient indicating that the stimuli 
were not of sufficient strength to ventilate the overload 

of co 2 . 

5. Alveolar 0 No overshoots were observed in relatively 

rapid on-transients. The half-times for the off-transients 
were longer than for the on-transients. There were over- 
shoots in the of f-trans ient s for the higher concentrations 


of CO^ input. 
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6. Transient half-times: On-transient half-times increased 

I 

O 

for f, V, and for increasing levels of CO^. On- 

transient half-times for P. decreased as the C0~ level 

°2 ■ 

was increased.. On-t rans'ient half-times for P 

1 co 2 

remained fairly constant for the various 
stimuli levels. In the off-transient V and V 

half-times decreased for increasing CO^ levels. Of f-transient 

P_ half-times increased for increasing C0 p levels . Once 
u 2 

again the half-times remained constant for increasing 

°°2 

stimuli CO^ levels.- Generally speaking, the off-transients' 
half-times were shorter for each variable than the corre- 
sponding on-transients except for P . The respiratory 

°2 

frequency seemed to have the most sluggish response 

of all the variables . 

O O 

Average response plots of , f, and V were given. The V 
versus time, t, plots for the various levels of C0 2 stimuli 
correspond to the results of similar simulations with 
Grodins * model although the exact inspired fractional 
concentrations of gases were not used. 

7. Subjective effects: Subjective effects included dyspnea, 

warmth, dizziness, headache, and visual disturbances all- 
with varying degrees of severity and at different levels 
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of stimuli. All subjects returned to a "normal" feeling 
after removal of the stimuli. 

The above summaries of the responses of the variables 
to 'the CO^ stimuli illustrate that much information can be ’ 
gained from studying such transient responses. These types 
of relationships are particularly important when subjects 
are exposed to rapid changes in the gaseous mixtures of the 
air. Certainly this type of experimental evidence is needed 
to establish 'a more physiologically sound model. 

F. A Model Evaluating CSF Subsystem Involvement in Ventilation 
*{• 

The H concentration is a variable that is well established 

as a contributor to the control of ventilation. Grodins* 

+ , 

model uses special weightings of the H concent-rations m 
the CSF compartment, the venous blood of the brain, and at 
the -carotid bodies' site to control ventilation.. With 
Grodins' model any combination of weightings for the 
CooWtt + \ or C rtr + \ ( t- T'_ 1 ) can be utilized. Thus the research 

L/bJ; ^ XI ) 3, \ il ) SLD 

of Milhorn et al . ( 6 )‘ which attempts to establish an appropriate 

detecting site within the brain structure can be used to 
justify the weightings m Grodins’ model. The primary purpose 
of this model was to establish a satisfactory transient 
response in ventilation which would apply if the CO^ detection 
site was assumed 'to he on either side of the central receptors. 

The controlled system aspect of the model was based upon 
the previous developments of Milhorn eL al . (l) and Grodins 
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et al. ( 9 ) vith the compartments being the brain,, bungs 1 , 
cereb'rospinal fluid,, tissue, and peripheral sensor (playing 
the role of the caroti'd bodies). The mathematical develop— 
merit is essentially the same as for Grodms 1 model except 
for the peripheral sensor. Since this is the most charac- 
teristic aspect of the model and the one which distinguishes 
the model from all others the equation which defines the CO^ 
tension at the proposed sensor site is repeated here. 


5 = P + ( P 

CS C0 2 b C0 2 . CSF C0 2 


where 



exp 



(15) 


B 


CO, 


= Deep brain tissue C0 ? tension, 


CSF 


CO, 


cs 


= Cerebrospinal fluid C0 p tension. 


= Average blood flow at sensor sit< 


S = CO^ dissociation curve's slope, 

D = Brain tissue CO^ diffusion coefficient, and 

X = Depth of sensor below medulla's surface. 

With this model ventilation is controlled by peripheral 
sensor H + concentration (C p0 ) and central sensor H + 


H + 


concentration (C cs ). In summary, C nR is given by 


H+ 


CS 


H+ 


) cs = exp (-2. 303 PH C0 ) 
H + 


where 


- n cs 


6.11 + lo-^’ 


C CS - K'O /Yf5o)p cs 
Cb C0„ b °2 G °C0. 


W co o /Y60) p cs 
2 


( 16 ) 


(17) 


CO. 


I 
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vit'h Of 


CO, 


and C 


CS 


CO, 


Solubility of C0 2 in brain compartment 
Concentration of COg at sensor site. (6) 


Pulmonary ventilation (v) is the sum of two signals. 

One signal involves the functional dependence upon C 

CS H + 

and the other signal -is a function of C—^ . The weighting 

b H+ 

is distributed such that for a given weighting of the central 
sensors the peripheral sensors’ signal has a weighting of 
100$ minus the central sensors’ weighting. Tidal volume 

. \ O 

t V T ) is then generated from a functional relationship of V. 

In a like manner dead space volume (V.^) is generated from a 
linear relationship with V . Respiratory frequency (f) is 

O 

the ratio V/V^, which is then multiplied by to yield minute 
dead space' (f\’^). The difference between minute dead space 

. O 

and pulmonary ventilation is the alveolar ventilation (V^J * 
Alveolar ventilation is then used as an input to the controlled 
system . 

The unique contribution of this research is associated 
with the determination of the sensor depth, X, which appears 
in Equation 15. A simulation was run for an inspired gaseous 

V 

concentration of 7$ COg with varying values for 'X. The range 

of sensor depth tested was o $ X S 1000P with the ventilation 

transient responding more rapidly for the sensor sites’ 

greater depths. A value of 280p was then used for further 

evaluation of the model. Utilizing a CO,, tension gradient 

this depth corresponds to a position which is 77$ of the 

distance between cerebjospin.il fluid P and deep brain 

°2 


CO, 


tissue P 
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Responses for three different controller equations 
were evaluated.. They were the linear controller equation of 
Milhorn (l) , 


O 





(18) 


another controller equation developed hy Grodins (9)> 

9 A - C CS1T +b 2 C a„ ’ 191 

n+ H + 

and the proposed controller equation of Milhorn et. al . (6), 


V = vf. (C_„ ) + (1-W) f (C ). ' (20) 

1 CS H+ 2 P ' S H+ 

Although no experimental data were shown the study indicates 
that the results of Equation 19 varied more from the normal 
expected response than the results of Equations 18 and 20.- 
When these three controller equations were compared for 
perfusion of the medulla surface with a mock cerebrospinal 
fluid COg tension of 71 mm Hg Equation 20 was indicated as 
giving the best representation of tidal volume (V^) response. 
Several other reported research efforts by Milhorn et al . (6) 

supported the simulation results for the CSF perfusion 
studies, the fact that CSF ventilatory control shouldn't 
he based entirely upon COg levels or H+ concentrations, and 
the fact that there is in all probability an optimum sensor 


depth . 
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■III. CONCLUSIONS 

A review of respiratory control system models developed 
"by Milhorn and hi3 associates has been presented. Chronologically 

I 

speaking, the models that were discussed encompass those 
developments associated with the steady-state models published 
In the roid~1960’s to the more recent investigation concerning 
a more detailed description of the CSF involvement with the 
control of respiration. The later developments also attempt 
to simulate transient responses as well as steady-state 
responses . 

The model utilized In the evaluation of the CSF sub- 
system’s function in the control of respiration most closely 
parallels the model developed by Grodins . (6,9). As was 

suggested in the Introduction one of the goals of this Study 
B'eport is to propose modifications for Grodins' model which 
stem from Milhorn' s work. Most of the suggested modifications 
involve some aspect of the controller equation. Any of the 
suggested modifications will not involve the exercise simu- 
lation subroutines that presently exist in Grodins ' modified 
program. The control of respiration and the system's variable 
responses that are related to exercise will need to be adjusted 
so that any new modifications are made compatible with 
existing exercise controller functions. 

One form of a controller equation that would be relatively 
easy to incorporate into Grodins' model is the one suggested 
by Milhorn with the special weightings indicated by Lambertsen. 
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(2) This equation is defined as 

V * 0.12 V (C PQ ) + 0 . UU V (C_ ) + 0.44, V (C ,) (21) 

PS H+ B H+ CSF H+ 

with all the variables defined in the same manner as with 

Grodins' model,. The neurological exercise contribution could 

be added to this equation. 

If Grodins' model was refined to simulate the actual 
gas exchange mechanism then the efforts of Milhorn and Pulley 
which describe the venous admixture phenomenon could be 

i 

implemented. With the existing Grodins' model no quantitative 
benefit can be obtained from the venous admixture description. 

If the ventilatory controller equation is considered to 
be a direct function of the 0 ^ and CC^ arterial gas tensions 
instead of the H+ concentration detecting systems then the 
additive and multiplicative relationship 


V = 2 (P - 37. 24) (1 + 13..6/(P - 25)) * 0 

C °2 ■ V, 


( 22 ) 


could be used. One must keep in mind that this equation is 

developed for P. £ 37.24 mm Hg. With the present 

00 

version of Grodins^ program the alveolar CO^ tension doesn't 
remain above 37.24 mm Hg when simulating variations in exercise 
under the environmental conditions of 5 Torr CO^ at 260 mm Hg. 
Consequently, the limitations of this equation will need to 
he redefined’. The variables involved here are easily accessible 
ip Grodins program. Comparison of 'the system's responses 
using a fox'd of this eouatJ on to the responses obtained - 
using the weightings of eompar tmonhai H+ concentrations could 
he beneficial in establishing appropriate detector sites. 
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Also, the relationships which were developed in conjunction 
with Equation 22 between dead space, respiratory frequency, 
alveolar ventilation and minute- ventilation (some of 'which 
are already us-ed in Grodins * program) could be implemented. 

In particular, a more, physiologically satisfying description 
of respiratory frequency could be obtained. 

A more complex modification would involve the sensor 
site location for the COg tension as described by Equation 15 
and repeated here. 


T rs = p b + “ p n ) ex P 

CS C0 2 B C0 a csf C0 2 b co 2 



(.2 3) 


Since the sensor depth (x) has been experimentally justified 
by Milhorn et al. (6) and all other terms are common to 
Grodins 1 program, the components of Grodins' controller 
equation which rely upon the H+ concentration of the venous 


blood of the brain and the CSF subsystem could be replaced by 


a form of Equation 23. This would be one component of the 
minute ventilation and would effectively combine the sensing 
mechanisms of the brain and CSF compartments. If it is 
more satisfying to have the controller equation be a weighted 
function of H+ concentrations the P can be converted into 

cs co 2 

a term corresponding to the H+ concentration at the effective 
sensor site. This relationship is also developed by Milhorn 
et al. (6) 


As previously implied, any modification to Grodins' 
program evolving from Milhorn' s research will need to be 
adjusted so that the exercise phenomenon can he simulated. 
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The latter two modifications that are suggested would be the 
most beneficial to pursue. They would add some detail and 
additional flexibility for adjusting responses which are 
missing in the present modified version- of Grodins 1 model. 
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6.3 Program Listing for the Modified Grodins' Respirator y Contr ol System 

Minor modifications are indicated in the right hand margins of subroutines RC13 and 
RC17. The addition to RC13 is to prevent the simulation from getting into a loop which 
was found to occur under certain conditions. The RC17 modification provides additional 
print routines as mentioned in the text. 
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